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TECHNICAL NOTE 3632 

CORRELATION, EVALUATION, AND EXTENSION OF LINEARIZED 
THEORIES FOR TIRE MOTION AND WHEEL SHIMMY 
By Robert F. Smiley 


SUMMARY 


An evaluation is made of the existing theories of linearized tire 
motion and -wheel shimmy. It is demonstrated that most of the previously 
published theories represent varying degrees of approximation to a sum- 
mary theory developed herein which is a minor modification of the basic 
theory of Von Schlippe and Dietrich. In most cases where strong differ- 
ences exist between the previously published theories and the summary 
theory, the previously published theories are shown to possess certain 
deficiencies . 

A series of systematic approximations to the summary theory is 
developed for the treatment of problems too simple to merit the use of 
the complete summary theory, and procedures are discussed for applying 
the summary theory and its systematic approximations to the shimmy of 
more complex landing-gear structures than have previously been considered. 

Comparisons of the existing experimental data with the predictions 
of the summary theory and the systematic approximations provide a fair 
substantiation of the more detailed approximate theories. However, some 
discrepancies exist which may be due to tire hysteresis effects or other 
unknown influences. Thus, further work may be needed to explain these 
discrepancies . 


INTRODUCTION 


In the ground maneuvering of aircraft equipped with swiveling landing 
gears there sometimes arises the problem of violent oscillations or shimmy 
of the landing gear which may lead to failure of the gear. In the past 
this problem has been handled largely by means of various measures based 
on practical experience. However, this empirical approach has not proved 
entirely adequate. Moreover, for radically different types of complex 
flexible landing gears it is highly doubtful whether any empirical approach 
based purely on past experience could always safely and optimumly take into 
account all of the possible conditions which a landing gear might be sub- 
jected to in actual operation. 
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A considerable amount of theoretical and experimental work on wheel 
shimmy has been done, mostly in the past 25 years. (Most of the existing 
papers on this subject are listed in ref. 1, which also presents a his- 
torical discussion of the development of the wheel-shimmy problem.) How- 
ever, most of these theoretical papers have not been .correlated with each . 
other or with the available experimental data, so that essentially there 
exists at present a large number of at least superficially different 
theories of wheel shimmy and a fair amount of experimental data which has 
not been correlated with many of these theories (refs. 2 to 23). 

The primary purpose of the present paper is to clear up this partial 
confusion of theories by demonstrating that most of the previously pub- 
lished theories represent various approximations to one basic general 
linearized theory derived herein and that most of the previously published 
linearized theories which do not represent approximations to this general 
theory possess certain undesirable characteristics. This basic general 
theory, which is henceforth called the summary theory, is derived in such 
a manner that it makes use of and is compatible with the soundest features 
of practically all the previously published theories, insofar as this ’is 
possible at present) however, in the main this summary theory is a minor 
modification of the theory proposed by Von Schlippe and Dietrich in refer- 
ences 3, 4, and 5- 

A second purpose of this paper is to develop a series of systematic 
approximations to the summary theory suitable for use in the treatment 
of problems too simple to merit the use of the complete summary theory 
and to examine both these systematic approximations and the previously 
published theories to determine how these theories are related to the 
summary theory and how the predictions of these theories agree with the 
available experimental data. 

A final purpose of this paper is to illustrate procedures for applying 
the summary theory and its approximations to complex types of flexible 
landing-gear structures. 

Although the primary purpose of this paper is concerned with the 
wheel-shimmy problem, most of the material presented is directly appli- 
cable to the more general problem of the motion of elastic tires under 
arbitrary rolling conditions. Thus this material is pertinent to' the 
study of veering-off or ground looping, ground handling, and catapulting 
stability of aircraft . 

The material in this paper is arranged as follows. First, a detailed 
statement of the problem is given, together with a detailed outline of the 
manner in which it is treated herein. Then, after a brief discussion of 
the restrictions on the analysis, a linearized derivation is made for the 
general Von Schlippe-Dietrich type of kinematic equations governing the 
motion of elastic tires rolling without skidding. This analysis proceeds 
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essentially in accordance with the theoretical analysis of Von Schlippe 
and Dietrich except that the present analysis considers the subject of 
tire tilt in slightly greater detail. 

Next, the primary forces and moments acting on a rolling tire are 
discussed and used to establish the equations of motion for arbitrary 
rolling conditions. Then a systematic procedure is developed for forming 
approximations to the summary theory. 

The previously published theories are listed, discussed, and com- 
pared with the summary theory and these systematic approximations. 

Finally, the application of the summary and approximate theories to 
several simplified landing gears is discussed. The first example is 
chosen primarily to demonstrate the correlation between theory and experi- 
ment, the second example to demonstrate the correlation between the sum- 
mary theory and its systematic approximations, and the remaining example 
to illustrate the application of the theory to complex problems. 

Some of the material presented herein was submitted to the University 
of Virginia in partial fulfillment of the requirement for a Master of 
Aeronautical Engineering degree. 


SYMBOLS 


a trail (perpendicular distance between ground-contact center 

point and swivel axis) 

, . T , sin k 

a' = a - £Lh — - — 

Ai,Ag,A3 coefficients defined by equations (115b) 

Bi,Bg coefficients defined by equations (ll5b) 

c lateral distance of center of pressure of vertical force 

from XZ -plane 

c-^ change in lateral distance of center of pressure of vertical 

force from XZ-plane per unit of 7 \q 

c y change in lateral distance of center of pressure of vertical 

force from XZ-plane per unit of 7 

c-j- distance from wheel center to center of gravity of swiveling 

parts of landing gear 
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e 2 

d l> d 2 

D 

D( ) 


D t ( ) 

E 0> E 1> 

E 

f 

Wl' 

F yh 


F • 
yi 


yn 

F ys 

F yA 
F 


yy 


F r 


g 

h 


distance from center of gravity of swiveling parts of 
landing gear to swivel axis 

hysteresis constants used in equation (l28c) 


tire parameter used by Bourcier de Carbon 


differential operator with respect to distance , 


d ( ) 


or v -1 D, ( ) 


differential operator with respect to time, 


dx 


d ( ) 

dt 


or v D ( 


coefficients of linear differential equations 
energy dissipated per cycle 
frequency, v/2jt 

coefficients of linear differential equations 


lateral force due to hysteresis effects 

lateral inertia force resulting from lateral deformation 
of tire 

net lateral tire force acting on wheel 

net lateral structural force acting on wheel 

lateral force on tire due to lateral distortion of tire 


lateral force on tire due to lateral tilt of tire 
vertical load on tire 

lateral force on swiveling parts of landing gear due to 
landing-gear strut 

linear damping constant (Damping moment = g 
half-length of tire-ground contact area 
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i = \JZ I 

l( ) imaginary part of ( ) 


I 


o 


J yt 

I yw 

•^i(f 


moment of inertia of the swiveling part of a landing gear 
about an axis parallel to the swivel axis and passing 
through the center of gravity of the swiveling part 

polar moment of inertia of wheel and tire about an axis 
perpendicular to the wheel axle 

polar moment of inertia of tire (excluding solid wheel parts) 

total polar moment of inertia of wheel and tire 

moment of inertia of the swiveling part of a landing gear 
about the swivel axis 


j 

J 

k 

k n 

AKa 


excess of number of zeros over number of poles 
parameter in stability-determination plots (appendix C) 
lateral spring constant of landing-gear strut 
parameter used in appendix A 
torsional stiffness of tire 

total effective change in tire torsional stiffness due to 
tire inertia effects 



effective change in torsional stiffness of tire due to 
lateral acceleration of tire 


K 7 

k a 

a^a 


change in tire torsional stiffness due to centrifugal forces 

lateral tire force due to tilt per radian of tilt angle 
lateral stiffness of tire 

total effective change in lateral stiffness of tire due to 
tire inertia effects 


AK- 


*i 


effective change in lateral stiffness of tire due to lateral 
acceleration of tire 
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X Q y ~L ^ b s > 

L 

m 


“t 

“w 

m l 

M„ 


M. 


xs 


M 


X0 


M, 


zh 


M zi 


M, 


zs 


M, 


za 


M, 


z/ 


M 


zA 


Mq 



change in lateral stiffness of tire due to centrifugal 
forces 

(nL + h)h n_1 

tire constants; = 

' n n! 

relaxation length 

mass of swiveling parts of landing gear 
mass of tire 

mass of wheel including tire 

mass of nonswiveling parts of landing gear 

constant friction-damping moment 

net structural tilting moment acting on wheel center 

gyroscopic moment due to swiveling 

twisting moment due to hysteresis effects 

inertia moment resulting from lateral deformation of tire 

net structural swiveling moment acting on wheel center 

torsional moment on tire due to twist of tire 

gyroscopic moment due to tilting 

gyroscopic moment due to lateral distortion of tire 

torsional moment on the swiveling parts of the landing gear 
due to landing-gear strut and damper unit 

damping moment about swivel axis 

cornering power (lateral tire force per radian of yaw angle 
during steady yawed rolling for yaw angle approaching 
zero) 
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n l 

P 

Pl>P 2 


parameter used in appendix A 

complex roots of characteristic equations 
functions defined in appendix D 


p-i oo^Ppoo^P -i -| jP -| p ^P pi ^>Pgp functions defined in and after equations (80) 


q n parameter used in appendix A 

r free tire radius 

r polar radius of gyration of tire 

O 

r^ vertical distance from wheel axle to ground 

r^ radius of cross section of tire torus 

R tire parameter used hy Bourcier de Carbon 

R( ) real part of ( ) 

s circumferential coordinate on tire (fig. 1) 

S wave length, 2ft jv 

S tire parameter used by Bourcier de Carbon 

t time 

At time lag due to tire hysteresis 

T Moreland's time-lag constant 

T tire parameter used by Bourcier de Carbon 

- * * functions of B^- correlating structural forces, moments, 
and deflections 

T time-lag constant for hysteresis moment 

(JO 

T-^ time-lag constant for hysteresis force 

u = c-,F (aZ-i cos k - a)sin k 
k 7 \ z v ' 
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= c^F (l-j_ cos k - a) sin k 
v rolling velocity 

w width of tire -ground contact area 

w^ density 

w K = (aK^ - aF z + c^,F z sin K)sin k 

x horizontal distance parallel to mean direction of rolling 

motion 

X,Y,Z space-fixed coordinate axes; the X-axis is horizontal and 

parallel to the mean direction of rolling motion, the 
Z-axis is vertical, and the Y-axis is perpendicular to the 
XZ-plane. The XY-plane is the ground plane. 

y lateral distance of tire equator from XZ-plane 

z vertical distance up from XY (ground) plane 


lateral deflection of tire equator from XZ-plane; subscript 0 
refers to the center of the ground-contact area, 1 to the 
foremost point of the ground-contact area, 2 to the rear- 
most point of the ground- contact area, i to equator points 
off the ground, and g to equator points on the ground 

a twist in tire, radians 

Ax half -width of twisting-moment— angular-deflection hysteresis 

loop of tire 

7 lateral wheel tilt, radians 

7^ lateral tire tilt resulting from lateral deformation, radians 

e pneumatic caster, K^R 

- 

T lo' 1 ll> T l3 

■ lateral deflection of center plane of wheel with respect to 

n a >V the XZ-plane; subscript 0 refers to the point corresponding 

J to the center of the ground-contact area, 1 to the point 
corresponding to the foremost point of the ground-contact 
area, 3 to the center point of the wheel, a to the point 
of attachment of the swiveling parts of the wheel to the 
swivel axis, i to wheel-plane points off the ground, and 
g to wheel-plane points on the ground 
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V 

4 Z 


inertia-force parameter (eq. ( 46 )) 
inertia-moment parameter (eq. ( 49 )) 
hysteresis-moment parameter (eqs. (62)) 
hysteresis-force parameter (eqs. (60)) 


^2’^i ,Ag 


A A r 


\ =( £ 


angle of rotation of wheel about the vertical Z-axis, radians 
inclination of swivel axis, radians (fig. 5) 


lateral distortion of tire equator with respect to the solid 
parts of the wheel; subscript 0 refers to the center of 
the ground contact area, 1 to the foremost point of the 
ground contact area, 2 to the rearmost point of the ground 
contact area, i to equator points off the ground, and g 
to equator points on the ground 

half -width of lateral-force— lateral -deflection hysteresis 
loop of tire 


circular frequency of shimmy motion, 


2rtf or v^v 


w 


path frequency of shimmy motion, 

tire tilt parameter (eq. (13)) 

spring constant for a linear restoring moment 


P„ = S aK y - aF z + ac ^ F z + c 

a = 1 + Slh - K 

l ir 


yF z sin 


K^sin k 


a l’ a 2 


T 

T, 


constants representing phase shift 

constant defined by equations (33) and (50) 
constant for gyroscopic moment 


NZ-j^T 

~ 
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angle of rotation of wheel about the swivel axis , radians 
to angular velocity of wheel about its axle ^ 


Subscripts : 
c critical 

max maximum 


STATEMENT OF THE PROBLEM AND GENERAL APPROACH 


The purpose of this section is to define specifically the problem 
considered in this paper and to clarify further the correlation between 
the various parts of the paper. 

The basic problem to be considered is the rolling motion and wheel 
shimmy of a rigid wheel equipped with an elastic tire, when the wheel 
is attached to some supporting structure such as a landing-gear strut. 

The motion of the rigid wheel can, of course, be completely described 
by six independent variables corresponding to the three degrees of free- 
dom in translation and rotation of the wheel. In addition to these 
six degrees of freedom, there exists a seventh degree of freedom which 
is associated with the distortion of the elastic tire or the track of 
the tire on the ground which results from the application of a given 
motion to the rigid wheel. Thus, in general, the motion of a rigid wheel 
with an elastic tire represents a system of motion involving seven varia- 
bles, and seven equations correlating these different variables are 
required to solve for the motion of a landing gear under arbitrary rolling 
conditions. Six of these equations will usually be the equations expres- 
sing the sum of the forces or moments acting along each of the three 
principal coordinate axes; the seventh relation will be an equation, 
usually a kinematic equation; which correlates the tire distortion with 
the other variables. 

The present paper is not concerned with all seven degrees of free- 
dom. Most of the paper is restricted to a consideration of cases of 
wheel motion in which the wheel is rolling at an approximately constant 
velocity v without braking, and consequently with constant angular 
velocity n>, and where no strong vertical oscillations are involved. 

Thus, for example, effects of acceleration or deceleration, which are 
known to have at least some influence on the rolling motion (see, for 
example, the experimental evidence of ref. 1J) are not considered. 
Similarly, fore and aft oscillations of the wheel are excluded. 
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When these three restrictions are applied* the seven-variable prob- 
lem of a rolling wheel becomes reduced to the consideration of a system 
involving the following four degrees of freedom: (l) swiveling of the 

wheel about a vertical axis through the wheel center point., designated 
by the symbol 0; (2) lateral tilting of the wheel with respect to a 
vertical plane parallel to the direction of undisturbed motion* desig- 
nated by the symbol 7; (3) lateral displacement of the wheel with 
respect to a space-fixed reference axis parallel to the direction of 
undisturbed motion* designated by the symbol tj with various subscripts; 
and (A) lateral displacement of the tire footprint on the ground (which 
is a measure of the tire distortion)* designated by the symbol y^. 

(These coordinates and their positive directions are illustrated in 
fig. 1.) 

In order to obtain four equations correlating these four variables 0 * 
7* T]* and y^* the following procedure is used: After some remarks on 

general restrictions* a kinematic relation between the four variables is 
derived in the section entitled "Kinematic Relations for the Rolling 
Tire." Next* the primary forces acting on the wheel from the ground* 
including wheel inertia forces* are discussed in the section entitled 
"Forces and Moments on the Wheel." By utilizing these ground forces and 
moments* the four basic equations of motion for the wheel* including the 
kinematic equation* are set down in the section entitled "Equations of 
Motion." 

For many applications these equations of motion in their most gen- 
eral form are relatively complicated and* although they are by no means 
insolvable* it is profitable to simplify the equations for those problems 
which do not require the detailed equations of the summary theory. There- 
fore* a number of systematic approximations to the summary theory are 
formulated in the section entitled "Systematic Approximations to the 
Summary Theory." A second reason for establishing these systematic 
approximations lies in the fact that they furnish a framework for com- 
paring the summary theory with the other existing theories of wheel 
motion* most of which are closely related to these systematic approxima- 
tions. Such a comparison of the summary theory and its systematic approx- 
imations with the existing theories of wheel motion is carried out in the 
section entitled "Classification and Evaluation of Existing Theories." 

In the last major section of this paper the summary theory and its 
systematic approximations are applied to three illustrative types of 
landing-gear configurations which are chosen either to illustrate agree- 
ment between theory and experiment or to illustrate methods for applying 
the theory to complex problems of wheel shimmy. 
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GENERAL RESTRICTIONS 


Before entering upon the detailed derivation of the equations of 
motion, some further restrictions on the analysis should be discussed. 
First of all, the present paper is limited almost exclusively to lin- 
earized theories. However, there is some question as to whether a lin- 
earized theory is sufficient to describe the important features of wheel 
shimmy. It appears at present that a linearized theory will provide at 
least a fair qualitative description of stability boundaries for shimmy 
and will indicate whether a given motion is stable or not. However, 
agreement between theory and experiment, presented in a subsequent sec- 
tion, is still not good enough quantitatively to warrant the conclusion 
that nonlinear effects can always be neglected or replaced by equivalent 
linear effects. 

Another limitation of the linearized theory is that it does not 
permit calculation of the maximum steady-state shimmy amplitude for those 
steady-state self -excited shimmy motions which sometimes occur on actual 
landing gears. 

Although the preceding considerations suggest that nonlinear effects 
in landing-gear motions may possibly be of importance for some practical 
problems, their consideration is beyond the scope of the present paper 
and henceforth only linearized theory is discussed. The only concession 
to nonlinearity is made in appendix A, which presents a conventional 
approximate method for converting a nonlinear shimmy damper to an equiva- 
lent linear damper. It should, however, be noted that some attention has 
been given to the development of nonlinear tire-motion theory in refer- 
ences 18 to 21. 

Another restriction arises in connection with the assumption adopted 
throughout this paper that the finite width of the tire need not be taken 
into account in developing a tire-motion theory for single tires of con- 
ventional cross section. This assumption appears at present to be at 
least partly justified on the basis of an experiment by Von Schlippe 
and Dietrich (ref. 3); on the other hand, since their investigation of 
this matter was extremely limited in scope, their experimental result 
may not be completely typical. Consequently, a more thorough evaluation 
of tire-width effects seems desirable. Some theoretical work on this sub- 
ject has been done by Von Schlippe and Dietrich (ref. 3) and later by 
Rotta (ref. 2), but the matter is beyond the scope of the present paper. 


KINEMATIC RELATIONS FOR THE ROLLING TIRE 


In this section the kinematic equations for the motion of a rolling 
tilted elastic tire without skidding are derived in accordance with the 
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theoretical analysis of Von Schlippe and Dietrich (ref. 3 or A) . This 
derivation differs only slightly from that analysis in that it omits 
some refinements of the theory which are necessary for very wide tires 
and it includes some influences of tilting of the tire in more detail. 
While the modifications that are made in regard to tilt may not neces- 
sarily he of practical importance in most cases, they may he of interest 
in a few problems. 

Specifically, the object of this section is to obtain a relation 
correlating the absolute lateral deflection of the center point of the 
tire ground- contact area Jq with the corresponding wheel coordinates 

of lateral deflection q (for example, t)q or tj^) , swivel angle 0 , 

and tilt 7 . (See fig. 1.) First, some geometric relations are set 
down and some background information regarding tire distortion is dis- 
cussed. Then this information is utilized to obtain a kinematic rela- 
tion between the lateral deflection of the tire center line or equator 
at the forward edge of the ground-contact area y^ and the coordi- 
nates rj, 7 , and 0. Next, a kinematic relation between the lateral 
deflections of the tire equator at the center and forward edge of the 
ground-contact area (designated yQ and y-^, respectively) is estab- 
lished. These two relations are combined to obtain a basic kinematic 
equation correlating y^ with tj, 7 , and 0 . 

The derivation of these kinematic relations is based upon the fol- 
lowing physical concept: As a tire moves forward, the tire material on 

the circumference just ahead of the ground-contact area is laid down or 
developed on the ground without skidding and becomes the new forward por- 
tion of the ground-contact area, so that the track of the tire is com- 
pletely determined by the lateral-distortion coordinate of the foremost 
ground-contact point y^ and the slope of the distorted center line or 

equator of the tire at that point. 


GEOMETRIC RELATIONS 


The primary geometric quantities involved in the problem of a 
rolling tire are shown in figure 1 , which gives an instantaneous view of 
a distorted tire with respect to an arbitrary space-fixed XYZ coordinate 
system, the X-axis being horizontal and parallel to the mean direction 
of wheel motion, the Z-axis being perpendicular to the ground, and the 
Y-axis being perpendicular to the X- and Z-axes. Parts (a) and (b) of 
this figure represent side and bottom views, respectively, of a rolling 
wheel that is swiveled and tilted. For the sake of clarity, part (c) of 
this figure, which shows an end view of the rolling tire, has been drawn 
to a different scale from part (b) and represents the unswiveled condition. 
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In discussing the geometric quantities, the following terminology and 
symbols are used: The wheel center plane is the plane of symmetry of 

the wheel perpendicular to the wheel axle. The tire center line or equa- 
tor comprises the tire points which on the undistorted tire are located 
at the intersection of the tire outer circumference with the wheel cen- 
ter plane; under the action of moments and lateral forces these tire 
points are deflected laterally by an amount A with respect to the wheel 
center plane. The symbol Aj_ designates the lateral deflection of tire 

equator points which are not in contact with the ground and A g desig- 
nates the lateral deflection of points which are in contact with the 
ground. The point at the center of the ground-contact area is designated 
by A 0 . 

The lateral distance of the wheel plane from an arbitrary space- 
fixed XZ-plane is designated by rn for points off the ground at a ver- 
tical height z and by q g for points on the ground. The lateral dis- 
tances of tire-equator points from this XZ-plane are similarly designated 
by and y g . The difference between y and tj is the lateral dis- 

tortion A of the tire, or 


Ai = yi - % (l) 

and 

Ag = Xg ~ \ (2) 

The tire contacts the ground in a finite area having a finite width 
and a length 2h. The width of this area is assumed to be negligibly 
small; that is, the ground-contact area is assumed to be reduced to a 
ground-contact line. The foremost ground-contact point (in the direction 
of motion) is designated by the subscript 1, the rearmost point by the 
subscript 2, and the center point by the subscript 0. Except for 
braking and accelerating effects, the center point 0 has approximately 
the same horizontal x-coordinate as the wheel axle. 

Distances about the tire equator or circumference are measured in 
terms of the circumferential coordinate s whose origin is taken at the 
point 0. 

The wheel is assumed to move at constant velocity v approximately 
in the direction of the X-axis. The -wheel ii laterally inclined with 
respect to the vertical Z-axis by the tilt angle 7 and is swiveled with 
respect to the XZ-plane by the swivel angle 0. Both tilt and swivel 
angles are assumed to be small; that is, cos 0 *=» cos 7 « 1, sin 0 « 0, 
and sin 7^7. 
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The center point of the wheel axle is located at a vertical dis- 
tance r^ from the XY (ground) plane , a lateral distance r 7 from 

the intersection of the wheel plane and the XY-plane, and a lateral 
distance T)* from the XZ-plane, where 

r \ 3 = t} q " (3) 


TIRE DISTORTION 


This section contains a short discussion of the features of tire 
distortion which are pertinent to the derivation of the basic kinematic 
relations of this paper. 

Experimental and theoretical considerations (for example, see 
refs. 3 and 2, respectively) indicate that, if the tire equator in the 
ground-contact region is subjected to arbitrary lateral distortion, the 
lateral distortion of the tire equator off the ground A^ tends to die 
out as an exponentially decaying function of the circumferential dis- 
placement s (for example, see fig. 2(a)). Thus, near tire point 1 
off the ground the tire distortion will tend to approach the pattern 
described by the equation 

s-h 

\ = v L (4) 

and a similar equation will apply near tire point 2. The exponential 
constant L is a tire characteristic having the dimension of length 
and is called the relaxation length. The relaxation length near point 2 
is not necessarily exactly the same as that near point 1; however, since 
the former relaxation length will not be used in this paper in any criti- 
cal calculations, this difference will not be taken into account. 

In regard to the accuracy of equation ( 4 ) very near point 1 , it 
should be emphasized that this exponential variation is only an expression 
of the equilibrium condition which the tire-equator distortion would 
reach in the absence of any restraints. However, it is obvious that con- 
ditions exist for which this distortion curve cannot be completely expo- 
nential in form. For example, for the case of pure lateral deflection of 
a stationary tire, the tire equator in the ground-contact zone is (neg- 
lecting skidding) a straight line parallel to the wheel center plane and 
extending from point 1 to point 2 (see solid lines in fig. 2(b)). Con- 
sequently, the existence of an exponential curve just to the right of 
point 1, and including point 1, would imply the existence of a sharp 
bend in the tire at point 1 such as is indicated in figure 2 (a) . Since 
a sharp bend is impossible because of finite tire stiffness, it follows 
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that, in general, on a stationary tire the exponential variation given 
by equation (4) cannot be valid close to point 1. However, experimental 
evidence indicates that beyond a short transition region ahead of point 1 
the tire-equator distortion curve does have an essentially exponential 
character (see solid lines in fig. 2(b)). As the wheel rolls ahead the 
nonexponential transition region of the tire equator is laid down or 
developed on the ground as it passes into the ground-contact zone, and 
the more nearly exponential part of the equator curve moves down toward 
the ground (see dashed lines in fig. 2(b)) and is eventually developed 
on the ground, so that after rolling a short distance from rest and 
during normal rolling conditions (fig. 2(c)) the tire-equator distortion 
at the front end of the tire can approach the assumed exponential varia- 
tion of equation (4). 

At the rear end of the tire the equator distortion curve during 
rolling does not so closely approximate an exponential variation, since 
at the rear end there is no process of laying down or development such 
as is responsible for the exponential variation at the front end. How- 
ever, since the rearward section of the tire equator is not used in any 
critical calculations in this paper, its equator curve is also, for sim- 
plicity, assumed to be exponential. 

If equation (A) is accepted as the basic equation for tire-equator 
lateral distortion near point 1 under rolling conditions, the total lat- 
eral displacement of the tire f!rom the XZ-plane in this region can, by 
use of equation (l), be written in the form 

_ s-h 

y ± = % + L (5) 

Substituting the geometric relation Hi = % - 7 Z ( see fig* 1) into 
equation ( 5 ) gives 

s-h 

= Hg - 7Z + \e L (6) 


KINEMATIC EQUATION 


By making use of the physical concepts discussed previously, together 
with equation (6), it is now possible to establish as follows the basic 
differential equation relating the tire deflection at the center of the 
ground-contact area y-Q with the wheel coordinates tj, 0, and 7. 

There is assumed to be perfect adhesion between tire and ground, 
that is, no skidding. As the tire rolls forward (arbitrarily swiveling. 
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tilting, and moving laterally) a distance dx, a new element of the tire 
of circumferential length ds above and in front' of point 1 is laid 
down or developed on the ground. This tire element, before being laid 
down on the ground, had the lateral-distortion variation given by equa- 
tion (6). This equation, after differentiation with respect to s, 
yields for a given instantaneous position of the tire the following rate 
of change of distortion: 


Hi = lk_ 7 iE.i Al e~ S£h 

ds ds ds L 


( 7 ) 


At point 1, where s = h and y^ = y-jj 


(Hi) . (Ik) . r (<k) . 1 A 

\ds / 1 \ds 1 \ds Li 


( 8 ) 


The term 


ds/- 


is simply the sine of the angle between the ground and 


the tire equator at point 1. (See fig. 1.) Just to the left of point 1 

the tire is flattened on the ground, or — = 0. If (— ) were not 

\ds/. 


ds 


zero, the tire would have to have a sharp bend at point 1. However, 
because of the finite bending stiffness of an actual tire, a sharp bend 

is impossible; thus = 0 and equation (8) reduces to 

(M) = - I A 

\ds ), Vis /j, L A 1 


( 9 ) 


Further, since = 0, s is a horizontal coordinate near point 1. 

The rate of change of wheel lateral displacement tj with respect to 

o 

the horizontal coordinate x at any given instant is just the swivel 
angle 0; hence 


(Hi 

Vds . 


= e -ih 


(10) 


If the tire is assumed to have no sharp bend at point 1, 


(ar), ■ (as 2 ), at thiB point - Then ' Blnoe (2 s ), 


is the slope of the 
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tire equator on the ground at point 1 and since no skidding is assumed 
to exist, this slope must coincide with the track of the rolling tire 

dy-i 

on the ground, which is — Thus, 


5i„ e .1 

dx L 


or, if differentiation with respect to x is designated by the operator 

D = — and the terms are rearranged, 
dx 


L Dy 1 = L9 - Ap (ll) 

(Alternate derivations of this equation are presented in refs. 3 and 4 .) 

A slightly more convenient form of equation (ll) is obtained by substitu- 
tion of the geometric relations A-^ = yp - T]p and T|p = T io + k® (see 

fig. l) to give 


(1 + L D)y = q 0 + (L + h)0 (12) 

Equation (12) is the basic equation for y 1 previously obtained by 

Von Schlippe and Dietrich (refs. 3 and 4 ). It should be noted, however, 
that no tilt terms appear in the equation. Although it is not known 
whether the effect of tilt on the validity of equation (12) is important, 
in view of the present lack of a reliable method for taking this tilt 
effect into account, the following argument is presented to afford at 
least a crude approach to the problem. 

Equation (8) contains the tilt term which was set equal to 

zero on the grounds that the factor is zero because of the finite 

bending stiffness of the tire. (See fig. 1.) On the other hand, if it 
is assumed that the bending stiffness of the tire is zero and if radial 

tire distortion is neglected, will be equal to h/r and the tilt 

term yh./x will enter into equations ( 9 ) to (12) . As a somewhat ques- 
tionable approximation, it will now be assumed that a term of this type, 
but smaller by the reduction factor £ < 1, should appear in the dif- 
ferential equations ( 9 ) to (12) . Equation (12) then becomes 

IhL 


(l + L D)y 1 = ti 0 + 0(L + h) 


r 


7 


( 15 ) 
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A similar tilt term was derived by Greidanus (ref. 7) on 'the basis 
of a slightly different argument. (Greidanus 1 term is discussed in a 
subsequent section of this paper.) However, apparently no other detailed 
tire-motion theory has included such a term. 

Equation (l3) is the fundamental kinematic relation for tire point 1. 
The kinematic relations for points 0 and 2 are determined by the condition 
of perfect adhesion between tire and ground. During the rolling process 
each tire circumferential element first contacts the ground at point 1, 
later proceeds to point 0 and then to point 2, after which it leaves the 
ground. Consequently, with perfect adhesion each tire element at point 2 
has the same lateral deflection that it had when it entered the contact 
zone at point 1 a distance 2h ago; that is, 

y 2 ( x ) = y x (x-2h) (1*0 

Similarly the kinematic relation for point 0 is 

y 0 U) = ypU-M 

or 

y x (x) = y Q (x+h) (15) 


Finally, by combining equations ( 13 ) j (l5)> and (5) the equation 
[l + L D]y 0 (x+h) = t} q (x) + [L + h]e(x) - ^ y(x) 

= Tj^(x) + [i. + hje(x) + jr^ - i^j 7 (x) ( 16 ) 

is obtained. This is the basic kinematic equation correlating the tire 
lateral deflection y^ under rolling conditions with the swivel angle 0, 

the wheel lateral displacement tJq or and the lateral tilt 7 for 

arbitrarily applied variations of 0, t|q or T)^, and 7. However, this 

transcendental form of the kinematic equation is not the most convenient 
form for some purposes in this paper. In particular it is expedient to 
remove the transcendental expression from equation (l6) by use of a 
series expansion. 
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SERIES EXPANSION OF KINEMATIC EQUATION 


The expression, y^x+h), after expansion in a Taylor series, gives 

2 

y Q (x+h) = y Q (x) + h Dy Q (x) + D 2 y Q (x) + . . . 


- y 0 ' x > + H £ D V X > 


si 111 


(17a) 


3 n 


n 

where the operator D represents — — » An alternative form of equa- 

dx 


tion (17a) which is useful later is 


y Q (x+h) = e^y^x) (l7b) 

since the infinite series in hD is the series expansion of the exponen- 
tial function. A third useful form of this equation is obtained by 
expressing equation (17b) in terms of a time derivative instead of 

the space derivative D. Since it is assumed throughout this paper that 
the rolling velocity v is constant, the correlation between these two 
derivatives is given by the equation 


D t ( ) » 12. = 2 12 . y 12 = v d( > 
dt dt dx dx 

and hence equation (l7b) can also be written in the form 


-1 


y 0 (x+h) = e hv y Q (x) 

Differentiation of equations (17) gives the result 
Dy 0 (x+h) = Dy 0 (x) + h D 2 y Q (x) + i h^y^x) + 


(17c) 


00 , n-1 

h 


n=l (n-1 ) ' 


7 D n y 0 ( x ) 


(18a) 
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or 


kD _1 hv "4) + 

Dy 0 (x+h) = De y Q (x) = v D t e 


y 0 ( x) 


(18b) 


Substitution of equations (17) and (18) into equation (l6) gives after 
rearrangement^ with y^(x) written simply as y^ and similar treatment 

of ^ 3 ’ 0 > and 


lo + l i e - 7 - 1 5 + V + ( r 3 - ^ 

= + D + *2 D 2 + . . .)y Q 

*( 1 + g l n D “) y o 

where 

. I « L + h 

l 2 = (2L + h) | 

l = (nL + h) — — 
n v ' nl 

or 

+ l i e - ^ 7 = ^3 + h Q + ( r 3 “ ^) 7 

= (1 + L Dje^ 

= (l + Lv -1 D t ^e hv ^ y Q (19b) 

Equations (l9a) and (19b) are alternative forms of the basic kinematic 
equation (l6) which are useful in subsequent sections of this paper. 
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This concludes the derivation and discussion of the basic kinematic 
equation correlating the lateral tire deflection Yq with the wheel 

coordinates 0, and 7 . Next, attention will be directed to the 

relationships existing between these coordinates and the forces and 
moments acting on the wheel. 


FORCES AND MOMENTS ON THE WHEEL 


In this section the primary forces and moments acting on a rolling 
wheel are discussed and, where possible, equations are set down for these 
quantities. These equations are then utilized in later sections, together 
with the preceding kinematic equation, to establish the equations of 
motion for a rolling wheel. 

The forces and moments considered fall into five general categories: 
elastic forces and moments due to tire distortion, gyroscopic moments, 
tire inertia forces and moments, hysteresis forces and moments, and 
structural forces and moments. 

Throughout this discussion, forces along the coordinate axes are 
considered positive if they tend to move the wheel in the positive direc- 
tions of the coordinate axesj moments about the coordinate axes X, Y, 
and Z or other parallel axes are considered positive if they tend to pro- 
duce wheel rotation from the positive Y-axis toward the positive Z-axis, 
from the positive Z-axis toward the positive X-axis, and from the posi- 
tive X-axis toward the positive Y-axis, respectively. 


ELASTIC FORCES AND MOMENTS DUE TO TIRE DISTORTION 
Lateral Elastic Force 


The lateral elasticity properties of a tire will be considered first. 
If a static untilted tire is laterally deflected at its base with respect 
to its rim by a lateral force Fy^, it produces an equal spring reaction 

force roughly proportional to the mean lateral distortion , or, 

inversely, a lateral tire distortion T y^^ creates a proportional 
ground force Fy^. If the lateral distortion of the center of the ground- 
contact line 7 \q is taken as the mean distortion, then the elastic ground 
force is 


F y* - K x\> - - %) - %(yo - 15 - r 3 r ) 


( 20 ) 
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where is the lateral spring constant or side stiffness of the tire. 

This relation is used by most investigators. However, in references 2 
to 5 a slightly different expression is used. In these references the 
mean lateral distortion of the tire is defined as the average of the 
distortions at the leading-edge and trailing-edge points of the ground- 
contact area (points 1 and 2) . The resulting equation for F ^ is 

F yX “ I K x(h + \>) < 21 > 

instead of equation (20). The true equation for Fy^ is probably more 

complicated than either of these two equations; however, since no plausi- 
ble means of obtaining a better equation is available, it appears advis- 
able to select one of the above equations for use in this paper. Equa- 
tion (21) may be slightly the better equation for a few special cases of 
wheel motion, but equation (20) is much simpler to work with, and in most 
cases of wheel motion it makes little difference which of the two equa- 
tions is used. Therefore, for the sake of simplicity equation (20) is 
adopted hereinafter as the basic equation for the lateral force on a 
wheel due to lateral deformation of the tire. 


Torsional Elastic Moment 

The torsional elasticity properties of a tire will be considered 
next. If a tire is twisted on the ground about a vertical axis through 
an angle a, there arises a restoring ground moment that is roughly 
linearly proportional to the twist : 

“za = V < 22 > 

The tire twist a is equal to the mean angle between the track of the 
tire on the ground and the wheel plane; that is, a = Dy - ean - 0. Taking 
the value of Dy mean as Dy Q gives 

a = Dy 0 - 0 (23) 

and thus 

“za - XaK - 6 ) - ^'Vo ' 6 ) (2k ~> 

Most investigators of tire motion use this relation. However, in refer- 
ences 2 to 5 the mean angle is taken equal to - Ag ^2h and thus 

the moment equation 



( 25 ) 
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is obtained, which leads to relatively more complicated equations of 
motion than does equation (2^). Since there is no strong reason for 
believing equation (25) to be a significant improvement over the simpler 
equation (2^), the latter is used in the analysis of this paper. 

Melzer (ref. 10) has used the less accurate relation that the moment 
due to tire twist is 


M m = -1^0 (26) 

which implies the relation 0 » Dy Q . (See eq. (24).) Since this rela- 
tion is not true in all practical cases, Melzer' s theory should be viewed 
with some caution. 


Tilt Elastic Force 

If a tire is tilted from the vertical Z-axis by an angle 7 without 
lateral distortion of the equator (Aq = 0 ), there arises a restoring 

ground lateral force that is approximately linearly proportional to the 
tilt angle (e.g., see ref. 2); 


F y7 = -K 7 7 (27) 

where is the constant of proportionality. Most authors (excepting 

Rotta in ref. 2) have not considered the effects of this force term 
although they have considered other effects of the same order of magnitude. 


Vertical-Load Center of Pressure 

Under some circumstances the vertical load F z influences the 

wheel motion. In order to consider this influence it is necessary to 
know the location of the center of pressure of this force. In the 
XZ-plane (fig. 1) this center of pressure lies approximately below the 
wheel axle in line with the point 0. In the YZ-plane the center of pres- 
sure is shifted laterally from the intersection of wheel plane and 
ground tjq as a result of lateral distortion Ag and tilt 7. As a 

first approximation, this shift may be taken as linearly dependent on 
Aq and 7 so that the lateral distance c of the center of pressure 

from the XZ-plane becomes 
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o - no + c a\> - V 
‘ c A y o + i 1 - 0 a)% - °r 7 

“ °A y O + ( x - c a)">3 + p 1 - °A) r 3 - c r] 7 (28) 

where and are constants. (The signs of the terms are chosen 

so that c^ and c^, are positive numbers.) 

GYROSCOPIC MOMENTS 

Gyroscopic Moment Due to Lateral Distortion of Tire 


The origin of gyroscopic moments on a rolling untilted wheel with 
lateral distortion of the tire at the ground (fig. 3) is considered next. 
While the solid rim and axle parts of the wheel are untilted, lateral 
deformation of the elastic tire causes the tire, on the average, to be 

Vi 

tilted with respect to the wheel center plane by an amount 7^ = T + 


where r is the tire radius and t 1 is a correction factor which indi- 
cates the effective fraction of the total tire mass that is tilted at this 
angle. Kantrowitz (ref. 8), apparently the only investigator who has con- 
sidered this at least theoretically significant factor, has 

as 1/2. This tilting action produces an angular velocity 
where D^. indicates differentiation with respect to time. 


suggested that 

D 7 _ D tVl 
U t?A r + ry 

This angular 


velocity, together with the rotational velocity of the tire o>, produces 
a gyroscopic moment about the Z-axis of magnitude 


M zA - -V" Va < 2 9) 

where 1^ is the moment of inertia of the tire (excluding the solid 

rim and axle) about the wheel axle. By using the relation D^ ( ) = v D( ), 
equation (29) can also be expressed in the form 

M zA - 7 » 7 A (5°> 

where the ratio v/co is, to a good enough approximation for this sec- 
ondary term, equal to the tire radius r. Then, substituting for 7^ 
and cu/v in equation (30) gives 
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M, 


zA 


T llytv 
r(r + r 3 ) 


DAr 


(31) 


For later convenience, the result can he expressed in several alternate 
abbreviated forms: 


M zA = " tv2dA 0 = -tv 2d (y 0 - tj 0 ) = -tv D t (y 0 - ri 3 - r j7 ) 

where 


T 


T l I yt 
r^r + r 


J) 


( 32 ) 


(33) 


Another method for deriving an expression for t is discussed in a 
subsequent section. 


Gyroscopic Moment Due to Tilting of Wheel 

If the entire wheel structure tilts at an angular velocity D^7, 
another gyroscopic moment arises of magnitude 


M sr ' - V D t r 


T~ D t r 


in addition to the term of equation (29) . Here I 


yw 


(34) 

is the total polar 


moment of inertia of the wheel (including the tire) about its axle. 


Gyroscopic Moment Due to Swiveling of Wheel 

If the wheel swivels at an angular velocity Dj.0, a tilting gyroscopic 
moment also arises of magnitude 


m x0 = -Iy^ D t 0 D t e (35) 

TIRE INERTIA FORCES AND MOMENTS 


This section is concerned with an examination of the influence of 
tire inertia forces and moments on a wheel rolling at high speeds. Two 
types of such inertia effects are evaluated now in separate sections: 
inertia forces and moments associated with lateral distortion and twisting 
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of the tire, and centrifugal forces and moments. Then the overall effects 
of these two types of inertia forces and moments are considered in another 
section. 


Inertia Forces and Moments Due to Lateral Distortion of Tire 

At high rolling or sh immy velocities, tire inertia forces and moments 
arise which are proportional to the relative accelerations of the different 
parts of the tire (including the previously discussed gyroscopic moment due 
to tire lateral distortion, discussed here from a slightly different point 
of view) . A rough estimate of these forces and moments can he made as fol- 
lows: One-third of the total mass of the tire m^. is assumed to be loca- 

ted on the periphery of the tire and to- be subjected to the same accel- 
erations with respect to the wheel hub as are tire particles on the equa- 
tor line, while the remaining tire mass is assumed to be substantially 
undisturbed. The "active" mass of the tire per unit circumferential 
length is then r. The lateral acceleration of tire particles on 

the right-hand side of the tire and off the ground in figure 1 (a) will 
be considered first. The lateral distortion of the tire in this region 
is given by equation ( 4 ). The lateral relative velocity of a tire par- 
ticle, obtained by differentiating this quantity with respect to time, 
is 

s-h 

Ah - (Ah - t D t s ) e L 

The quantity D^s, which represents the peripheral velocity of tire 

particles with respect to the wheel axle-, is approximately equal to the 
negative of the rolling velocity v, so that the velocity expression 
becomes 


Ah = ( D th + 1 b) e " l 

Differentiation of this result to give the relative acceleration of the 
tire particles yields the result 

/ s-h 

A% - (vt\ + §7 Ah + ^ h) e " L 

The corresponding inertia force AF for this part of the tire is obtained 
by integrating the product of this acceleration and the active mass per 

nn P S55 *r o 

unit length to obtain the force term / ®t ^i ^ s * Evaluation 

6 jtr J s=h 
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of this integral, after replacing the upper limit by infinity for simpli- 
fication of the result (which introduces no significant error because of 
the rapidly decaying exponential function in , yields 

^--g( D A + T D th + ph) (36) 

The corresponding inertia moment 2M is given by the expression 

r sin cp ds 

where r sin cp is the moment arm (see fig. 1(a)), s is related to cp 
by the relation s - h = r(cp - cp-^), and cp^ = sin ^ p. Therefore, the 
moment integral may be written in terms of cp in the form 

* ■ -§?£ r sl " ‘‘’ p 2 ’' 1 + ¥ + ? Ai ) e L r dtp 


^Js=h 


Evaluation of this integral, after replacing the upper limit by infinity 
(which introduces no appreciable error), yields the expression 



In a similar manner, for tire particles off the ground on the left- 
hand side of the tire in figure l(a), the following expressions are 
obtained for the inertia force and moment : 


2 M = 



V 2 



D t A 2 



(38) 

(39) 


In these two expressions it is assumed, for reasons previously discussed, 
that the relaxation length L is the same for both sides of the tire. 
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In obtaining the inertia forces and moments for tire particles in 
the ground-contact area, it is recognized that in practically all cases 
where inertia forces are important the ground-contact line is almost a 
straight line, so that the lateral distortion for tire particles in this 
region can be expressed fairly well by the equation 

Ag = Aq see 

The corresponding velocity and acceleration are 


D t A g = D t A 0 + s D t a - va 


D t% = D t% + s D t 2a " 27 D t a 
The total inertia force for this region is then 


AF 


ul os-h - 
= / D t Ag ds 


6jt 


r J s=-h 


'g(°tA 0 -2vD t 


a) 


and the inertia moment is 


AM =- 


“t 

n s=h 
1 

6itr J 

s=-h 

mth? 

D + ^a 

SWr 

t 


s d£ 


g 


(^ 0 ) 


(41) 


The total inertia force Fyp, obtained by summing the force terms 

in equations (36), (38), and (40) , can be stated conveniently in terms 
of Aq and a by using the relations Ap + A2 = 2Aq and Ap - A2 = 2ha, 

which are valid for a substantially straight ground-contact line. The 
result is 






where = L + h. Similarly, the total inertia moment is 


( 42 ) 
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The significance of these inertia expressions will now be partly- 
evaluated by considering the inertia force for sinusoidal oscillations; 
that is, where Aq = '^Q max s i n ^ and -therefore D + ^Aq = -v 2 Aq, so 
that equation (42) may be restated as 



__^t_(y£ 

3itr \L 



*0 




(An equation similar to equation (44) has been derived by Marstrand in 
reference 20. Marstrand 1 s equation, however, is based on a cruder repre 
sentation of the shape of the lateral distortion of the tire.) 


In order to interpret the significance of the inertia force it is 
noted that the tire force quantity which is of importance for the sub- 
sequent analysis is the net tire force F-^ acting on the wheel, which 

is equal to the sum of the ground force Fy^ and the inertia force F^ 


Fyn - Fy^ + Fy^ (45) 

For a static tire, F 1 ^ was set equal to K^Aq (eq. (20)). In the 

dynamic case the relation between ground force and lateral distortion 
of the tire may be modified by the inertia effect. As a first approxi- 
mation, it will be assumed that the modification of the ground force is 
proportional to the inertia force, or 


F yA “ K A A 0 - Vyi ^ 

where tj is a number whose absolute value will be less than unity if 

y 

the modification of the ground force due to the inertia force is less 
than the inertia force itself. 

After combining equations (44), (45), and (46), the following 
equation for the net tire force Fy n is obtained: 
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F yn - F yx + F yi - 




m. 


3rtr \ 4 



(^7) 


From the form of equation. (47) it can be seen that, insofar as the ratio 
of net tire force to lateral deformation is concerned, the effect of the 
inertia force can be considered equivalent to a change in tire lateral 
stiffness AK^_ equal to 



m 


Similarly, from an exa m ination of the terms containing a in the inertia 


moment equation (43), it can be 
inertia moment is to change the 
AKq^, which is defined by 


AKa.= 


concluded that part of the effect of this 
tire torsional stiffness by an amount 


i 1 - \) 

tmt. 

rh | 

[h + n 

k -$) 

3* 


j2 + r 2 


- Lv 


hV 


3r 


m 


where T)„ is a number representing the torsional stiffness similar to 

'Z 

T}y for the lateral stiffness. The remaining inertia-moment term in 

equation (43), which is proportional to is simply the previously 

discussed gyroscopic moment due to lateral tire distortion. By comparing 
this term with equation (32) it is seen that the coefficient t may be 
expressed by the equation 

ayfc + Ly/i-ig) 

3 *(l 2 + r 2 ) 

Equation (50) gives approximately the same result as 

Kantrowitz' assumption that The discussion 

range in which these stiffness changes are important 
after the effects of centrifugal forces have been considered. 


(50) 

equation (33) with 
of the velocity 
is postponed until 


Effects of Centrifugal Forces 

Another inertia effect that may become significant at high speeds 
is produced by the centrifugal forces acting on the individual mass ele- 
ments of the tire. These centrifugal forces appear to increase the tire 
stiffness, as will be demonstrated by a crude analysis which gives a 
qualitative idea of this effect but which should not be regarded as pos- 
sessing any strong quantitative merit. 
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For the purpose of this estimate, one-half the mass of the tire is 
assumed to be concentrated in the side walls and the other half is assumed 
to be concentrated on the periphery. 

If the tire lateral and torsional stiffnesses K A and are 

assumed to be directly proportional to the tension in the side walls of 
the tire, there will be two sources of tire stiffness: inflation pres- 

sure, which produces a side-wall tension approximately equal to wp per 
unit circumferential distance (where w is the tire width), and centri- 

^ corre- 
sponding to the peripheral tire mass ^ mj-. Thus the lateral stiffness 
of the tire may be expressed in the form 


fugal force, which produces the side-wall tension 


( 




or, equivalently, as 


K A <x ^itr^wp + rn^v 2 

K A = K A (l + -^1— \ 

A Astatic 

mtv 2 K A 


= k a 


static 


static 


4jcr 2 wp 


It is evident from this equation that centrifugal force increases the 
tire lateral stiffness by an amount AK A .: 

tJ 


m tv 2 K A . 
AK A . = Static 

J AjtrS/p 

and the torsional stiffness by an amount AK^.; 

J 


(51) 


,2 V 


mj.v -astatic 
AKq,. = - p — 

^ 4rt r^wp 


(52) 


Significance of Tire Inertia Effects With Respect 
to Tire Stiffness 

The significance of the two just discussed tire inertia effects on 
the tire stiffness will now be considered. 
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For the lateral stiffness of a tire, the effective change AK^ from 
the static value of is obtained by adding the two increments given 

by equations (48) and ( 51 ) . The resulting effective overall change in 
tire lateral stiffness as a function of rolling speed and shimmy fre- 
quency is 


ak a = 


l 1 ~ v) m t z iv 2 i 1 - yH' 


mj-v 2 ^ 


static 


3itr 


3jtrL 


4jtr^wp 


(53) 


The first term, which involves the shimmy frequency, appears to be small 
enough in comparison with so that it can probably be neglected for 


most practical conditions. The last two terms have opposite signs if 
t] < 1 and thus may represent two partly counterbalancing effects. The 

second term arises from the previous considerations of the lateral accel- 
eration of tire particles and tends to reduce the effective lateral stiff- 
ness of the tire with increasing rolling velocity if tj < 1. The last 

y 

term arises from the previous considerations of centrifugal forces and 
tends to increase the lateral stiffness. These last two terms indicate 
that at high rolling speeds, if r|y < 1, the tire stiffness may either 

drastically decrease or drastically increase, depending on which of the 
two terms is larger. However, both terms happen to be of the same order 
of magnitude and the derivations of both terms are based on concepts too 
crude to justify conclusions regarding which term is larger. Thus, the 
only conclusion that can be drawn is that at sufficiently high rolling 
speeds drastic changes in tire lateral stiffness may occur. Whether the 
stiffness increases or decreases can probably be settled only by 
experiment . 


In order to give some quantitative measure of the velocity at which 
these inertia effects become significant, some calculations were made to 
determine the velocity at which the magnitude of the second term in equa- 
tion (53) becomes equal to K^. By making use of the static tire data in 

reference 24 for several modern aircraft tires and assuming that Tjy = 0, 

it was found that this velocity averaged approximately 
400 /r fps « 270 sir mph, where r is expressed in feet . Similar esti- 
mates for the velocity at which the third term in equation (53) becomes 
equal to K\ yielded approximately this same velocity. Moreover, since 
this velocity is rather high compared with normal present-day landing 
speeds, the inertia effects on tire lateral stiffness considered here 
can probably usually be neglected. 


For the torsional stiffness of a tire, the overall effective change 
in torsional stiffness due to tire inertia and centrifugal forces 

is obtained by adding the two increments given by equations (49) and ( 52 ). 
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The result is 





static 


4rtr^wp 


( 5 ^) 


This equation is parallel to equation (53) for the lateral stiffness, so 
that statements made previously concerning the lateral stiffness apply 
here also. 


Other Inertia Effects 

The preceding discussion suggests that the effects of tire inertia 
are to change tire stiffness at high speeds and to introduce a gyroscopic 
moment. However, it should he recognized that other inertia effects will 
come into play, probably at velocities close to those at which the pre- 
viously mentioned inertia effects arise. For example, the basic kinematic 
equation depends on the assumption of an exponentially distorted tire- 
equator line corresponding to a definite "static" relaxation length. 

This assumption is valid (if it is valid at all) only when the elastic 
forces in the tire predominate over the inertia forces. Where inertia 
forces are strong in comparison with elastic forces, it is at least 
doubtful whether the relaxation length remains constant. 

Although there are undoubtedly other effects of tire inertia in 
addition to the ones discussed here, it appears probable that the impor- 
tance of many tire inertia effects can be assessed by means of the fol- 
lowing summary statement: The major effects of tire inertia on the rol- 

ling motion appear to come into play at a velocity of an order of magnitude 
of 400 \fr fps « 270 >/F mph where r is expressed in feet. For con- 
siderably smaller velocities, most inertia effects can probably be safely 
neglected) for velocities of this order of magnitude or higher, many of 
the basic assumptions of this paper, and of most other papers on this 
subject, may be subject to considerable error. 
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HYSTERESIS FORCES AND MOMENTS 


In addition to the forces and moments previously discussed, certain 
damping forces and moments arise as a consequence of the sometimes con- 
siderable hysteresis losses which arise in the distortion of elastic 
tires . Apparently the only significant attempt to deal with this hyster- 
esis problem is reported by Von Schlippe and Dietrich in reference 5 
This reference provides some valuable insight into the fundamental mech- 
anism of the hysteresis process and presents an equation for the hysteresis 
moment acting about the swivel axis of a shimmying wheel. However, even 
though use of this hysteresis-moment equation leads to good agreement 
between theoretical and experimental stability boundaries for a limited 
amount of experimental data (as is shown subsequently in the present 
paper), some parts of the analysis seem so unrealistic that it is ques- 
tionable whether much confidence can be placed in the final results of 
reference 5 - Apparently, the only other significant contribution to the 
hysteresis problem is provided by the analysis of Moreland in refer- 
ences 11 and 12. In these references tire hysteresis forces as such are 
not considered, but the idea is introduced that a tire possesses a char- 
acteristic time-lag constant. In a subsequent section of the present 
paper it is shown that this time-lag constant may be, at least in part, 
a consequence of hysteresis effects. However, the interpretation of 
Moreland's time-lag constant as a hysteresis effect presents some ques- 
tionable features that are also discussed subsequently. 

No completely satisfactory solution of the hysteresis problem has 
been found yet. However, the following crude analysis of this problem 
offers another point of view with a few qualitative merits not possessed 
by the two previous analyses. 

Consider the case in which a standing tire is subjected to a period- 
ical lateral deformation Aq of the form 


An = An sin vt 
u max 

Under these conditions the lateral ground force on the tire is 

experimentally observed to vary with time in the manner indicated in 
sketch 1 and the corresponding variation of lateral ground force F 


^Although Von Schlippe and Dietrich considered hysteresis effects 
in an earlier paper (ref. it-), this earlier analysis leads to some con- 
clusions which are not in agreement with the results of the later, more 
detailed analysis of reference 5« 
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with lateral tire distortion Aq > shown in sketch 2, appears in the form 
of a typical hysteresis loop. As can be seen from sketch 1, the lateral 




Sketch 1 


Sketch 2 


tire deformation lags behind the lateral ground force F y ^ by a 

time increment At, where At is approximately equal to the ratio of 
AAq to the m axi mum slope of the curve of Aq plotted against time 

(which is v for the assumed variation of Aq): 

\ max / 


At 



A 0 v 
u max 


( 55 ) 


K \ 

-AWWV- 


-H 


Sketch 3 


As a first approximation for quantitatively including this time-lag 
concept in the present analysis, a tire is 
assumed to behave somewhat like a combina- 
tion of a linear spring and a damper unit 
such as is indicated in sketch 3> where the 
spring constant corresponds to the previously 
discussed tire lateral stiffness and a^_ 

is the coefficient of an equivalent linear 
damper. Inertia forces are neglected for the 

present argument. The differential equation for this system is 




•A 


K A A 0 + a l D t\) - F yA 


(56) 


and its solution for the case of An = An sin yt gives a hysteresis 

u v max 

loop of the form indicated in sketch 2 where the time lag At becomes 


tan' 


a-, v 


At = 


-1 fl_ 

Ka 


(57) 


After equating equations (55) and (57) > a q can be expressed by the 
relation a-, = so that equation (56) can then be 

JL V 


V 
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written in the form 

F yX * K ?A> + p yh = h.{\> + Ta »t*o) (58) 

where Fy^, the lateral force resulting from hysteresis effects, is 


F yh = K A T A »t\) = K A T A V DA 0 


( 59 ) 


and where 


t a 


tall ( M °/ A 0 

V 


\ 


= tan 


^®max 



J 


(60) 


With the same type of reasoning the hysteresis twisting moment is 
given by the equation 


^zh — ^a^a, ®t a — ^a?-a v ■^ a 
(compare with eq. (59) )> where 

„ tan (*V°hax ) 

a v v 

► 

■n = tan 

a max 

It is now seen that the determination of the hysteresis force and 
moment from equations (59) and (6l), respectively, depends on the deter- 
mination of two quantities T^ and T a which have the dimension of 

time and which will be called time-lag constants. The quantity T^, 
in particular, can be considered closely analogous to Moreland's time- 
lag constant. In the present case equations have been derived for the 
time-lag constants as functions of the three variables t^, t^, and v. 

(See eqs. (60) and (62).) However, it should be clearly recognized that 
equations (60) and (62) are based in part on arguments valid only for a 
standing tire. These arguments may no longer be valid for a rolling 
tire, and even if the idea of a time-lag constant is still valid, it is 
likely that the time-lag constants will not be adequately predicted by 
equations (60) and (62), particularly if the quantities t) and t] z 

*7 

are evaluated from static hysteresis loops. Moreover, it is basically 
unsound to assume that the hysteresis force is dependent only on the 
tire lateral distortion Aq and is independent of the tire twist a| 


(61) 


(62) 
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actually the hysteresis force (and the hysteresis moment) will in gen- 
eral depend in a complex manner on both Aq and a, and even for a first 

approximation the interaction of these two variables cannot necessarily 
be neglected. Thus it appears that the preceding hysteresis equations 
are based on rather speculative and perhaps unsound assumptions, at least 
from a quantitative point of view, and for this reason these equations 
will not be incorporated into most of the derivations in subsequent parts 
of this paper. On the other hand, the preceding derivation may be suf- 
ficiently plausible to give some idea of the order of magnitude of hys- 
teresis effects, particularly since Moreland has indicated in reference 12 
that his experimental data (mostly unpublished) demonstrates the exist- 
ence of a time-lag effect in tire motion; consequently, in a few parts of 
this paper some mention will be made of the consequences of introducing 
the hysteresis force and moment terms that have just been derived into a 
wheel-shimmy analysis. 


STRUCTURAL FORCES AND MOMENTS 


The preceding discussion covers the major ground forces and moments 
and the gyroscopic moments acting on the wheel. In addition, forces and 
moments are exerted on the wheel by the supporting structure. These will 
be designated as Fy g for the net structural force parallel to the Y-axis, 

M xs for the net structural lateral tilting moment, and M zg for the net 

structural swiveling moment. These forces and moments include shimmy 
damper moments, spring restoring moments, inertia forces in a landing- 
gear structure (exclusive of the wheel inertia force), and spring forces 
arising from the flexibility of a landing-gear strut or of the fuselage 
of an airplane. In general, most of these forces and moments can probably 
be considered approximately linear except shimmy damper moments; however, 
even these moments can be replaced as a first approximation by equivalent 
linear damping moments. (See, for example, appendix A.) 

Within the scope of a linear theory, these structural forces and 
moments will depend in a linear manner on the wheel-center coordinates 

0 , and 7 according to expressions of the type 


Fys = Tl(D t ), 3 + T 2 (D t )e + Tj(D t ) 7 

(63) 

M xs - + I 5 (Dt) 6 + *6( D th 

(64) 

M z S = UWt + T e( D t) e + T 9 ( D th 

(65) 
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where the T' s are functions of the differential operator D^., sometimes 
called transfer functions, whose specific forms will depend on the type 
of landing gear in question. 

This concludes the discussion of the forces and moments acting on a 
rolling wheel. Now these quantities will be utilized to set up the basic 
equations of motion for a rolling elastic wheel. 


EQUATIONS OF MOTION 


DERIVATION OF THE EQUATIONS OF MOTION 


In this section the linearized equations of motion for a rolling 
elastic wheel are set down with the aid of the equations from the pre- 
ceding sections. 

The sum of the lateral forces acting on the wheel parallel to the 
Y-axis is set equal to the inertia reaction to give (see eqs. (20) 
and (27)) 

F ys + K *(yo " ^3 “ r 3?) “ V = “v D t\ (66) 

or, rearranging, 

F ys + Vo “ ( K A + “V " K r 3 + S ) 7 = ° (6T) 

The first term in equation (66) is the structural force, the second 
term is the net force on the wheel resulting from tire elastic and inertia 
forces flQ = _ + AKk where AK^ is given by equation (53) \ 

the third term is the lateral ground force resulting from tire tilt^ and 
m^ is the mass of the wheel (including the tire) . For reasons previously 
discussed, hysteresis forces and moments are not included either in this 
equation or in the following equations. 


Setting the sum of the lateral tilting moments about the wheel center 
equal to the inertia reaction gives (see eqs. (20), (27) > (28), and (35)) 


V + F z Wo + ( F “ ^3 + i 1 - c *) r - " c 


V x 3 " 7 


7 - H 3 


IWo - ^3 - r 3 7 ) - v] r 3 “ “"F - D t 0 = J xw W 


( 68 ) 
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or 

M X s + ( K A r 3 + F zC A )y 0 “ [V 3 + F z c^Jt) 5 - D t © - 

K A r 3 2 + K 7 r 3 “ F z r 3 C 1 “ C A) + F z c 7 + J xw D t 2 7 = 0 (69) 

The first term in equation (68) is the structural moment, the second 
term is the moment resulting from the vertical ground load, the third 
term is the moment of the ground forces resulting from tire lateral dis- 
tortion and tilt, the fourth term is the gyroscopic moment resulting from 
the swiveling motion of the wheel, and 1^ is the moment of inertia of 
the wheel about an axis through its center parallel to the X-axis. 

Setting the sum of the swiveling moments about the wheel center equal 
to zero yields the equation (see eqs. ( 24 ), (32), and ( 34 )) 

M zs + Ka^^yo " 9 ) “ tv D t (y 0 - ^ - ^ 7 ) - D t 7 = D^e (70) 
or 

Mss + (KaY' 1 - rv)D t y 0 + tv Vi} - ("t + I„ - -rr^jy D t? = 0 

( 71 ) 

The first term in equation ( 70 ) is the structural moment, the second term 
is the ..net moment resulting from tire elastic and inertia forces exclusive 
of the gyroscopic moment due to tire lateral distortion f K_, = + AK 

v y S"CEL"tXC CX 

where AKq, is given by equation ( 54 )), the third term is the gyroscopic 
moment resulting from tire lateral distortion, and the fourth term is the 
gyroscopic moment resulting from wheel lateral tilt. 

Equations (67), (69), (71), and (l6) or (19), together with the three 
auxiliary equations (63) to (65) , are the basic equations of motion for an 
elastic wheel. If the T-functions in equations (63) to (65) are known for 
a particular landing gear, these equations can be solved simultaneously to 
determine the rolling behavior of the gear. 

Next the question arises as to the most profitable method of solution 
of these equations for practical landing-gear problems. Essentially, the 
choice is between exact and approximate solution of the equations. In the 
past, exact solutions (omitting some of the less important terms previously 
mentioned) have been made only for the simplest case of a rigid swiveling 
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landing gear attached to a rigid fuselage (refs. .2, 4, and 5)« Although 
the exact solution of these equations for more complex problems does not 
appear to present any insurmountable difficulties, relatively complex 
transcendental equations may be involved, so that it is worthwhile to 
examine the possibility of finding simpler systematic approximations to 
the general equations. 

A second reason for investigating systematic approximations to the 
summary theory arises in connection with the correlation of the summary 
theory with the other existing theories. Superficially, in its present 
form, the summary theory does not closely resemble most of the other 
existing theories. However, the approximations that are presented sub- 
sequently make the correlations between the different theories fairly 
easy to see. 

Subsequent sections of this paper will be concerned with the problem 
of establishing a series of systematic approximations to the general equa- 
tions and the correlation of these approximations with the other existing 
theories of wheel motion. However, before proceeding with these two mat- 
ters it is convenient to digress slightly to consider the exact solution 
of the general equations for the case of steady yawed rolling, in order 
to establish several relations which will be useful in later sections. 


EQUATIONS FOR STEADY YAWED ROLLING 


For an untilted wheel which rolls at constant velocity at a constant 
small swivel or yaw angle, yg(x + h) = yg(x) = Constant, 0 = Constant, 

and Tjj = 7 = 0 , so that equations ( 2 ) (with y^ for (l 6 ), ( 67 ), 

and ( 71 ) reduce, respectively, to the relations 


II 

£ 

( 72 ) 

y 0 = (L + h)0 = 

(73) 

F ys + ^0 = 0 

(74) 

M zs - K o9 - 0 

(73) 


By combination of equations (J2) and (73) the tire lateral distortion is 
found to be 
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* 0 = IjQ (76) 

By combination of equations (73) and (74) the lateral force on the wheel 
is found to be 


F. 


ys 


= -w 


The quantity which represents the lateral force per unit yaw angle, 

is an important tire characteristic called the cornering power or lateral 
guiding characteristic of the tire. Later in this paper it is found con- 
venient to represent this quantity by a single symbol N, where 


N = l ± K^ (77) 

Another property of the steady yawed rolling condition that is of 
some interest is the distance of the center of pressure of the lateral 
force behind the center of the tire, which is sometimes called the pneu- 
matic caster e = -M zs /Fy S . This quantity, according to equations (73) 

to ( 77 ) > is equal to 


- ^ zs _ 

= "Fy S - N 


SYSTEMATIC APPROXIMATIONS TO THE SUMMARY THEORY 


(78) 


In this section the possibilities for simplifying the preceding 
equations of motion are discussed, and a series of systematic approxima- 
tions to the general equations of the summary theory is set down. 

All but one of the equations of motion (eqs. (l 6 ) or (19), (63) 
to ( 65 ), ( 67 ), ( 69 ), and ( 71 )) are usually simple linear equations and 
present no great difficulties. The exception is the kinematic equation, 
which was originally transcendental in form (eq. (l 6 )) and was later 
expressed as an infinite series of linear terms (eq. ( 19 a)). The most 
promising way to simplify the kinematic equation appears to be to assume 
that the series expansion in equation ( 19 a) is a rapidly convergent 
series in which all terms above a certain value of n can be neglected. 
The rapidity of convergence of the series and its significance cannot be 
fully determined without a knowledge of the particular landing-gear con- 
figuration considered. However, some insight into this question can be 
obtained by considering the ease of purely sinusoidal oscillations of 
iv-jX 

the form yQ = e , where the quantity is the path frequency. 

Substitution of this expression into the infinite series in yQ in equa- 
tion ( 19 a) yields 



MCA TN 3632 


^3 


i 1 + ]>_ l n D n )y 0 = (pi„ + i P2oo)y 0 

\ n=l ' 


(79) 


where 


■\ 


p i» - 1 - Vi + Vi 


p 2» - 1 1 V 1 - Vi + b v i 


(80a) 


Another form for the p's can he obtained by substituting the rela« 

into equation (l6) . The result is 


iv n x 

tion yQ = e 11 


= cos v-jh - Lvj sin v^h 
p = sin v^h + Lv-^ cos v^h 


(80b) 


The rate of convergence of the p series of equations (80a) can be 
tested for any given frequency by substituting numerical values of L, 
h, and v^_ into equations (80a) and (80b) and comparing the individual 

terns. A typical comparison is shown in figure 4 for the conditions 
L = 0.8r, and h = 0.5r. The abscissa of this plot represents the oscil- 
lation's wave length S = 2it jv-^ and the ordinate represents the p func- 
tions. The label p ip means that this curve represents the sum of the 
first two terms in the p^ series, and the other labels are analogous. 

(The approximation symbols will be explained later.) From this figure 
it is seen that the series converge very rapidly. From a purely quali- 
tative point of view the figure seems to indicate that, in dealing with 
shimmy wave lengths greater than approximately 4 tire radii, two terms 
in each series are sufficient to represent fairly well the exact varia- 
tions, for wave lengths greater than approximately 6 radii one term in 
the P2 series and two in the p^ series are sufficient, and for wave 

lengths greater than about 20 radii one term in each series is sufficient. 
(The wave lengths cited here represent only order of magnitude and are 
not necessarily quantitatively significant.) To correlate these observa- 
tions with the conditions of wave length likely to be encountered in prac- 
tice it can be stated that the experimental data of Von Schlippe and 
Dietrich (ref. k or 5) and Kantrowitz (ref. 8), which are probably fairly 
typical in this respect, demonstrate wave lengths which are about 4 radii 
long at zero rolling velocity and which increase with increasing rolling 
velocity. Thus it appears possible that the use of only a few terms in 
the series expansion may lead to a reasonable prediction of shimmy char- 
acteristics for practical operating conditions. 
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With the preceding considerations in mind the following approxima- 
tions to the general wheel-motion equations were established. 

APPROXIMATION A 

As a first approximation to the general kinematic equation (19a) all 
terns for n > 3 will be neglected. This gives the approximate differ- 
ential equation 

y 0 + z i % + l 2 + b = ^0 + h Q ~ t 1 7 

* + V + ( r 5 - (81) 

This equation, together with all the general force and moment equations 
previously discussed, is referred to hereinafter as approximation A. 

APPROXIMATION B 

A second, less exact approximation for equation ( 19a) is obtained 
by setting l n = 0 for n > 2. Thus 

y 0 + l 1 Dy 0 + l 2 D^y Q = tIq + IjQ - y 

= n 3 + l ± Q + (r 3 - (82) 

This equation is referred to as approximation B. 

APPROXIMATION Cl 

Another, cruder approximation for the general differential equa- 
tion (19a) is obtained by neglecting all terms in the series for n > 1. 
This gives the differential equation 

y 0 + l x Dy 0 = T] 0 + ^ 7 

= q 3 + l ± Q + ^r 3 - ^7 (83) 

which is referred to as approximation Cl. 
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APPROXIMATION C2 


As a slight simplification of approximation Cl, the relatively 
unimportant, or at least questionable, term involving | may be omitted 
from equation ( 83 ). This gives the differential equation 


y 0 + 1 D y 0 ■ 10 + 1 i 6 - 13 + ! i 9 + 


(»o 


which is referred to as approximation C2. 

With the aid of equations ( 2 ) and (23), equation (84) can be written 
in the more easily interpreted form 


Aq — — ^ 2^* 

or, by using in addition equations (20), (22), and (77 ) > as 


(85) 


_ _R_ a _ F yA _ ^za 

'° " K x ‘ K x 



( 86 ) 


Thus, in this approximation the lateral distortion of the tire is directly 
proportional to the angular distortion. 


The physical meaning of this approximation can be obtained by con- 
sidering that equation (86) can also be arrived at by letting the ground- 
contact half-length h approach zero in the general differential equa- 
tion (19a) (as was mentioned by Rotta in ref. 2), since all terms in the 
series for n > 1 and the tilt term are multiplied by h. Then equa- 
tion (l9a) (with 1=0) becomes 


y Q + L Dy 0 = q 0 + L0 

or, by using in addition equations (2) and ( 23 ), 


A = -La 

Also, equation (77) "the yawed rolling becomes 

N = K a L 

and the combination of equations (88) and ( 89 ) gives 


(87) 

( 88 ) 
(89) 


( 90 ) 
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Equation (90) is the same as equation (86) for any given combination of 
N and Thus, when written in the form of equation (86), approxima- 

tion C2 formally corresponds to the assumption of h = 0. 

Reliable qualitative results should be expected from approximation C2 
only when the neglected quantity h is small with respect to the char- 
acteristic length of the rolling motion in question (for example, the wave 
length s of a sinusoidal oscillation). Fortunately, this condition is 
at least sometimes satisfied for practical rolling conditions. 


APPROXIMATION D1 


Before considering the next approximation it should be remembered 
that all of the terms neglected in the preceding approximations were 
multiplied by the tire ground-contact half-length h; thus these approxi- 
mations implied the assumption of progressively smaller ground-contact 
length or progressively larger wave length. In order to simplify the 
equations further, it is necessary to make some assumptions about the 
other tire properties. Three such assumptions will now be made to sim- 
plify further the equations of approximation C2. For the first approxi- 
mation, to be called approximation Dl, the simplification 

l ± = 0 ( 91 ) 

is adopted. Then it follows from equation ( 85 ) that, for finite a, 

A 0 = 0 (92) 

which is the basic equation for this approximation. Thus for this 
approximation the tire is free to twist but not to deflect laterally. 
Infinite lateral stiffness is, therefore, also implied: 

k a = 00 (93) 

For the simplest form of wheel shimmy, due to tire elasticity rather 
than structural elasticity (to be considered subsequently), equation ( 92 ) 
does not provide accurate information. For wheel shimmy due largely to 
structural elasticity rather than tire elasticity, this approximation may 
be of some value; actually most existing theories corresponding to this 
approximation have been developed for the primary purpose of considering 
the influence of structural elasticity on wheel shimmy. 
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APPROXIMATION D2 


As a second simplification of approximation C2, the assumption 

= « (9*0 

can be adopted. The corresponding theory is designated as approximation D2. 
From equation ( 85 ) it is evident that this approximation implies that, for 
finite Aq, 

a = 0 (95) 


which in turn implies that 


N = 00 

= eN = 00 (96) 

Thus for approximation D2 the tire is considered to be torsionally rigid 
but laterally flexible. 


APPROXIMATION D3 


A third simplification of approximation C2 can be obtained by keeping 
the quantity l -j finite but considering the tire to have both infinite 
lateral stiffness and infinite torsional stiffness, or 


K A = = N = 00 (97) 

This approximation, which is designated as approximation D3, thus repre- 
sents the case of a rigid tire and consequently also implies that 
a = Aq = 0. (Formally, approximation D3'can also be interpreted as the 
limiting subcase of approximation D1 where N = °° or as the limiting 
subcase of approximation D2 where K A = However, it should not be 
concluded that approximation D 3 is necessarily inferior to these other 
two approximations . ) 

A choice of seven simplified approximations based on the summary 
theory is now available. The problem remains of determining which, if 
any, of these approximations is the simplest one which can be used for 
any particular tire-motion problem. While it is not yet possible to 
give a completely satisfactory answer to this question, some insight 
into the answer can be gained by comparing the various approximations 
with the previously published tire-motion theories, which are at least 
partly successful and most of which are closely related to these 
approximat ions . 
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CLASSIFICATION AND EVALUATION OF EXISTING THEORIES 


In this section the major previously published theories of wheel 
motion are briefly reviewed, evaluated, and, wherever possible, corre- 
lated with the preceding summary theory of this paper and its approxi- 
mations. Each of the major existing theories is first considered indi- 
vidually and then an abbreviated overall summary classification is 
presented in tabular form. 


INDIVIDUAL REVIEW AND EVALUATION OF EXISTING THEORIES 
Von Schlippe and Dietrich 


The tire-motion theory of Von Schlippe and Dietrich (refs. 3 to 5)/ 
of course, corresponds directly to the s umma ry theory of this paper, 
since the summary theory was taken from their theory with only minor 
modifications. These modifications include a more detailed consideration 
of some of the influences of lateral tilt and of tire inertia forces and 
moments. It should be noted, however, that the Von Schlippe-Dietrich 
theory is more advanced than the summary theory of this paper in that it 
partly takes into account the width of , the ground-contact area. However, 
as was previously noted, this factor is probably not of great practical 
importance . 


Rotta 

Rotta's tire-motion theory (ref. 2) corresponds to the summary theory 
of this paper because it also is based on the Von Schlippe-Dietrich theory. 
Rotta's theory represents a slight extension of the last theory to take 
into account more adequately most of the effects of tire tilt and the 
width of the gr ound - c ont a ct area. No inertia forces due to tire lateral 
distortion or centrifugal forces are discussed. 


Bourcier de Carbon Advanced Theory 

Bourcier de Carbon (ref. 6) has developed two closely related theories 
of tire motion which are similar to approximations B and C2. The first of 
these will be referred to as the Bourcier de Carbon advanced theory and the 
second as the Bourcier de Carbon elementary theory. 

4 

Bourcier de Carbon's advanced theory uses five basic tire properties 
which are correlated with those of the present paper through the following 
relations: 
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^9 


D = - 
N 

(98a) 

T = — 

(98b) 

k a 

S = -i- 

(98C) 


N 

(98d) 

2(L + h) 

h K a Vf 21 + h > 

(98e) 


Equations ( 98 ) were obtained by comparing this theory with the corre- 
sponding approximation B. The symbols of Bourcier de Carbon are over- 
scored and do not necessarily bear any relation to other symbols in this 
paper designated by the same letters. Although the symbols D, T, S, 
and ~e bear a simple relation to the symbols used in the derivations of 
the present paper, the symbol R bears a more complicated relation which 
is worth some detailed consideration. 

Bourcier de Carbon defined the tire property R as follows: If an 

untilted wheel is rolled forward while exposed to a constant turning 
moment about a vertical axis and with no side force, it will move in a 
circular path of a definite radius; R is defined as the reciprocal of 
the product of the turning moment and the path radius. Unfortunately, 
however, this constant -moment circle-rolling experiment is not easily 
performed. Therefore, equation (98e), which expresses R in terms of 
the more easily measured fundamental quantities L, h, and K^, is of 

importance for the use of the Bourcier de Carbon advaneed theory. 

In treating the subject of tilt, Bourcier de Carbon omits many of 
the details considered in this paper. For example, he implicitly assumes 
that Ky = c-^ = Cy = i = 0 and that the inclination angle k is small 
(taking cos k « l) . However, these omitted tilt terms may be as impor- 
tant as the terms considered (as will be shown later); therefore, Bourcier 
de Carbon's considerations of tilt are- incomplete. 

It should be noted that in reference 6 certain misconceptions occur 
in the parts of the paper that deal with comparisons between theory and 
experiment. In particular, some of the experimental data quoted by 
Bourcier de Carbon from reference 3 of the present paper appears to be 
either misquoted or misinterpreted. Consequently, Bourcier de Carbon's 
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conclusion that the experimental data of reference 3 provide a remarkable 
check of his theory is not completely justified; actually these experi- 
mental data provide only a fair indirect check of the theory. 


Greidanus 


Another theory similar to approximation B, except for the influence 
of tilt, is that of Greidanus (ref. 7)* Greidanus considers the influ- 
ence of tilt in much greater detail than does Bourcier de Carbon; how- 
ever, he also fails to consider the force term proportional to K^,, and 
thus his results also do not fully describe the influence of tilt. 

In addition, Greidanus' kinematic equation differs from equation (82) 
for approximation B in that he has introduced a slightly different term 
associated with tilting of the tire. In the present terminology Greidanus 1 
equation reads 

y 0 +l ! D y 0 + l 2 = TJ 0 + V - l 2r (99) 


The difference between the two equations lies in the coefficient of 7. 
For approximation B (eq. (82)) the coefficient is 


|Lh 

r 


( 100 ) 


and for Greidanus' equation (after substituting for 

( L * l) h 

r 


1 2 from eq. (19 a)), 
( 101 ) 


If £ is set equal to ^ j^L, the two coefficients are identical; 

thus Greidanus' kinematic equation can be considered to be a particular 
case of the corresponding equation of approximation B. 

No subsequent detailed discussion of Greidanus 1 theory is included 
in this paper because no complete translation of reference 7 is available. 


Bourcier de Carbon Elementary Theory 

Bourcier de Carbon's elementary theory corresponds to approximation C2 
of this paper except for the minor shortcomings which were discussed in con- 
nection with the Bourcier de Carbon advanced theory. The only difference 
in Bourcier de Carbon's two theories is that the coefficient R is taken 
as infinity in the elementary theory but i£ finite (see eq. (98 e)) for the 
advanced theory. The infinite value for R corresponds to the assump- 
tion l 2 = 0, which was previously made in passing from approximation B 
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to approximation C2 (compare eqs. (82) and (84)). The physical signifi- 
cance of R = 00 is obvious from equation ( 98 ). It means that h = 0. 


Melzer 

The Melzer theory of tire motion (ref. 10) is also similar to 
approximation C2 except for details of the tilting process. Otherwise 
Melzer *s kinematic equation is identical with the kinematic equation of 
approximation C2 and of Bourcier de Carbon's elementary theory. However, 
Melzer' s theory differs in that it treats the moment due to tire twist 
as proportional to the swivel angle -0 rather than to the tire twist 
angle DyQ - 0. This assumption would appear justified only if DyQ « 0, 

which is not true in general. It is interesting to note that for the 
simplest case of wheel sh immy (see section entitled "Application to Wheel- 
Shimmy Problems — Case I") the Melzer approximation leads to one of the 
same stability boundaries and to the same limiting high-speed shimmy fre- 
quency as the more nearly correct approximation that includes the term in 
DyQ. This restricted agreement, however, hardly justifies the use of 

Melzer 's approximation, since predictions made by means of the two approxi- 
mations differ with respect to divergence of the shimmy oscillations and 
with respect to another stability boundary. Moreover, for simple problems 
the Melzer approximation is not significantly easier to solve than the 
more nearly correct form including the DyQ term. 


Moreland Advanced Theory 

Moreland has proposed three versions of a tire-motion theory in 
references 11 and 12. The most advanced of these versions is governed 
by the equation 

a + TDta=-^=-M (102) 


or 


TZ-jV D 2 y Q + l ± Dy 0 + y 0 = 4o + TZ 1 V D9 + l l Q ( 10 3) 

where T is a time-lag constant. This theory corresponds to a general- 
ization of approximation C2 (with pneumatic caster neglected, that is, 
e = 0), since for T = 0 equation (102) is identical with the basic 
equation for approximation C2. However, for T 5 / 0 this theory is not 
directly compatible with the summary theory and its approximations. 
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Moreland uses the following reasoning to establish this equation: 
First, it is known that in steady yawed rolling a yaw angle a is devel- 
oped as a consequence of the application of a lateral force Fy A according 

to the relation 



(104) 


which is the basic equation for approximation C2. However, for the 
dynamic rolling case this equilibrium yaw angle obviously cannot be estab- 
lished immediately upon application of a given side forcej rather, a 
finite amount of time will be required for the equilibrium yaw angle to 
develop. Moreland has attempted to take this finite time lag into account 
by modifying equation (104) to the new form of equation (102). In the 
latter equation the constant T is a measure of the time lag of the yaw 
angle behind the applied force Fy^. 

This time-lag tern introduced by Moreland does not correspond exactly 
to any of the terms in the summary theory, and to this extent Moreland's 
advanced theory is apparently incompatible with the summary theory. How- 
ever, a partial reconciliation of the two theories can be obtained through 
the following considerations of hysteresis effects as applied to approxi- 
mation C2: According to equation ( 58 ) the lateral ground force, if tilt 

and inertia forces are neglected, is given by the equation 

F iv 

\> + t a D t*o - (“5) 

and the kinematic equation for approximation C2 is (eq. ( 85 )) 


A 0 — ^ 

Combining equations (105) and ( 85 ) to eliminate Aq and substituting 
K A li = N yields the equation 


a + T- 


A 


D t a = 


_yA 

N 


( 106 ) 


Equation (106) is formally identical with Moreland's basic equation (102) 
if the hysteresis time constant T A is considered equivalent to Moreland's 

time constant T. The important points to be noted here are: (l) according 

to both views, the tire twist a lags behind the applied lateral force Fy A , 

and (2) Moreland adopts the lateral-force equation Fy A = K a Aq which implies 
that the lateral force and lateral deformation are in phase (although lateral 
force and twist are not in phase), whereas according to approximation C2 
(see eq. ( 105 )) the lateral force lags behind the lateral deformation as a 
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consequence of the kinematic relation A = -Z^a. In regard to the second 

point, since there are apparently no pertinent experimental data available, 
it is not possible to conclude which point of view, if either, is corrects 

Why Moreland’s time- constant concept has not been incorporated 
directly into the derivations of the summary theory should be clear from 
the preceding discussion: It is not certain whether Moreland's time- 

constant terms are really independent of the terms already contained in 
the summary theory or whether they are, rather, another way of looking 
at some terms which are already included in the summary theory. More 
specifically, Moreland's analysis does not include inertia forces and 
moments due to tire lateral distortion, hysteresis forces and moments, 
or the higher l n terms (i 2 , Z^, . . . ), and for certain conditions 

any of these factors could be interpreted as a time-constant effect. In 
view of these factors and in view of the lack of pertinent experimental 
data, a completely satisfactory evaluation of the relative values of the 
summary theory and Moreland's advanced theory cannot be made in the 
present paper. 


Moreland Intermediate Theory 

As a simpler approximation for his advanced theory, Moreland has 
implied ,(ref. 11) that the influence of the time-lag term in the basic 
equation for his advanced theory (eq. (102)) can be approximated for the 
usual range of shimmy frequencies by using the simpler kinematic equation 

kOIja = -A 0 (107) 

Inasmuch as approximation C2 has the kinematic equation ( 85 ) 

1 2-CX = — Aq 

and approximation D2 has the kinematic equation (95 ) > which could be 
written in the form 


“ a -Aq 

it follows from a comparison of these last three equations that Moreland's 
intermediate theory falls somewhere between approximations C2 and D2. 

Since Moreland has not offered any concrete justification for this 
approximation, further detailed discussion does not appear warranted 
here . 
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Moreland Elementary Theory 

Moreland's most elementary theory corresponds directly to approxi- 
mation D3, the case of a completely rigid tire, except that it, like 
Moreland's other two theories, does not take into account the pneumatic 
caster (e = 0) , 


Temple Elementary Theory 

Temple has proposed an elementary theory for the motion of tires 
which is identical with approximation D1 (ref. 13)* Temple has chosen 
the most general form of this approximation in that he has considered 
both the tire torsional stiffness ■ (indirectly interpreted as an 

increase in trail) and the cornering power N. 

This theory was developed before experimental evidence pointing to 
the need for more detailed considerations of tire lateral stiffness was 
available. Subsequently, Temple has indicated a need for more refined 
considerations of the tire (ref. 25) and has developed independently a 
theory (unpublished, but partly described in ref. 21) similar to the 
theory of Von Schlippe and Dietrich. 


Maier 

Maier (ref. 14) has proposed a simplified theory similar to approxi- 
mation Dl, with the difference that he makes the added assumption that 
the tire torsional stiffness is zero. This theory, too, was devel- 

oped before there existed much experimental evidence pointing to the need 
for more refined considerations for shimmy behavior. 


Taylor 

Taylor (ref. 15 ), in a brief paper, suggested another tire motion 
theory which corresponds to approximation D2 except that details of the 
tilt process are omitted. 


Kantrowitz and Wylie 

The preceding theories for tire motion, which include most of the 
known theories, may all be considered as closely related to the summary 
theory of this paper. However, two other well known theories, one by 
Kantrowitz (ref. 8) and one by Wylie (ref. 9 ), apparently cannot be 
derived from the summary theory and thus cannot be accurately classified 
here with respect to the other theories. They possess some of the merits 
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of approximation B but in other respects are less adequate than approxima- 
tion C2. To point out the deficiencies of these' two theories it is suffi- 
cient to consider two simple cases of tire motion as follows: 


The first case to be considered is the steady straight-line motion 
of a nonswiveling, untilted, rolling wheel which is not yawed with respect 
to its direction of motion (a = 0) , which 
is inclined by an angle 9 (equal to the 
swivel angle) to the reference X-axis, 
and which has no lateral forces or moments 
acting on the wheel. (See sketch.) Obvi- 
ously, for this case there will be no lat- 
eral distortion of the tire, or 



A 0 = ° 


On the other hand, Kantrowitz' basic equation, which is 


A 0 + L DAq = 10 - l 2 D8 (108) 

gives for this steady unyawed case ^with DAo = D0 =0^ 

A 0 = L0 

This equation is obviously incorrect, since it implies that the lateral 
distortion of a straight rolling wheel, which actually must be zero, 
depends on the choice of the coordinate axes to which 0 is referred. 

Only for the special case of a wheel running along the reference axis 
(that is, for 0 = 0) is Kantrowitz' theory correct in this respect, and 
in an actual shimmy problem this condition is possible only for the case 
of zero trailj hence Kantrowitz' theory cannot necessarily be expected 
to give reliable results for trails different from zero. Thus Kantrowitz* 
theory is of at least doubtful value for practical shimmy calculations. 

In order to evaluate the Wylie theory, consider the case of steady 
untilted yawed rolling of a wheel moving parallel to the X-axis. (See 
sketch.) Obviously, the lateral distor- 
tion of the tire 7 \q will depend only 
on the swivel angle 9 (9 = a) and not 

at all on the absolute lateral displace- 
ment of the wheel tjq. On the other hand 

the basic equation of Wylie, which in the 
present terminology is 



y 0 + L Dy 0 = L0 - l 2 D0 


(109) 
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gives for this steady case (where DyQ = D0 =o) the relation yQ = IB 
or, by using equation (2), 


A 0 - L8 - tj 0 

This equation states the obviously incorrect conclusion that the tire 
distortion is dependent on tjq or, in other words, that it depends on 

the choice of the coordinate axes. Thus, only for the special case 

t) = 0 is Wylie's theory plausible in this respect, and tjq = 0 implies 

that the reference axis must pass through the path of the wheel. Since 
this condition is satisfied in an actual shimmy motion only for the 
special case of zero trail, Wylie's theory, like Kantrowitz', can be 
fully valid only for zero trail and, consequently, this theory is also 
of doubtful value for practical shimmy calculations. 

It might be noted that the preceding difficulty concerning Wylie's 
theory could be removed by the logical procedure of adding the term t)q 
to the right-hand side of Wylie's equation (109) to give the new basic 
equation 

y Q + L Dy 0 = IB - Z 2 D0 + ^ (110) 

Another questionable point in these two theories is that the 
Kantrowitz theory, as previously noted, predicts that the lateral dis- 
tortion in yawed rolling parallel to the X-axis (see previous sketch) is 

*0 = L0 

and so does the Wylie theory if the reference axis is chosen to give 
t)q *z 0. On the other hand, the summary theory has for this yawed case 

the relation 


A 0 = (L + h)0 

(See eq. ( 85 ).) The difference arises from the fact that Kantrowitz 
and Wylie did not explicitly consider the ground-contact length 2h in 
their derivation. 


Other Theories 

In addition to the just discussed theoretical papers dealing par- 
ticularly with the subject of landing-gear shimmy, a number of relevant 
papers exist which are either largely of historical interest, which deal 
particularly with automobile shimmy problems, which deal only briefly 
with landing -gear shimmy problems, which deal with other tire-motion 
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problems such as yawed rolling and veering-off or ground looping, or which 
deal with the determination of tire stiffness parameters. Although these 
papers are of some interest, they do not appear to contain any important 
contributions which are not contained in the theories just reviewed. The 
reader is referred to reference 1 for a substantially complete listing 
and brief discussion of most of the papers in this class. 

Of particular historical interest among the investigations not con- 
sidered here in detail are the work of Broulhiet (ref. 23) and the work ~ 
of Fromm (discussed in ref. 22). These two investigators independently 
were apparently the first to recognize the importance of lateral distor- 
tion and cornering power of tires in the wheel shimmy problem. Taking 
these factors into account, both authors developed tire-motion theories 
whose kinematic relations correspond to that of approximation C2 of the 
present paper. 


TABULAR CLASSIFICATION OF EXISTING THEORIES 


In order to permit easier visualization and comparison of the merits 
of the theories discussed, the major assumptions of the various theories 
of tire motion are collected in table I. This table lists the nature 
of the assumptions made in regard to the primary tire parameters N, 

K^, 1^, e, and Z n for each of the theories discussed. 

APPLICATION TO WHEEL-SHIMMY PROBLEMS 


In the preceding sections of this paper a set of basic differential 
equations for the motion of an elastic wheel has been derived and compared 
with the corresponding equations of most of the previously existing theo- 
ries. These comparisons have indicated that, from a mostly qualitative 
point of view, the summary theory of this paper and the systematic approx- 
imations to it incorporate the major merits of the existing theories of 
tire motion and avoid some of their disadvantages. However, it still 
remains to investigate how best to apply the theory to specific landing- 
gear problems, to investigate the question of the absolute or quantitative 
accuracy of the summary theory and of the other theories, and, if the 
s ummar y theory be found satisfactory, to establish the simplest systematic 
approximation to it which will give reliable information regarding any 
particular problem in tire motion. The best way to accomplish these vari- 
ous aims appears to be through a discussion of the shimmy of several par- 
ticular landing-gear configurations. In this section three particular 
landing-gear configurations are discussed which range in complexity from 
the simplest case of a rigid swiveling landing gear to the most general 
case of a gear of arbitrary complexity. 
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DESCRIPTION OF PARTICULAR CASES. CONSIDERED 


The first landing-gear configuration considered, which is designated 
case I, is illustrated in figure 5« It consists of a rigid landing gear 
whose only degree of freedom other than tire distortion is rotation of 
the wheel about an inclined swivel axis , which may be opposed by a linear 
spring or damper. This particular configuration is chosen because most 
of the existing experimental data have been obtained for such a configu- 
ration. Thus, this configuration makes it possible to discuss and evalu- 
ate the summary theory, its systematic approximations, and the existing 
theories with respect to agreement with experiment in regard to the var- 
ious important characteristics of a shimmy motion, such as stability 
boundaries, shimmy frequency, and divergence. 

The second landing-gear configuration, case II, is an untilted 
landing gear possessing two degrees of freedom aside from tire distortion. 
This landing-gear configuration, which is illustrated in figure 6, con- 
sists of a wheel free to swivel but not to tilt, which turns about a 
rigid vertical swivel axis attached by a spring k to the supporting 
structure. (This spring is an idealized representation of the lateral 
flexibility of an actual landing-gear strut.) This configuration is dis- 
cussed for two purposes: (l) to illustrate the effect of structural 

elasticity on wheel shimmy behavior and (2) to provide an example which 
is better suited than case I for bringing out the relative merits of 
several of the systematic approximation theories for a case involving 
structural flexibility. 

The third landing-gear configuration considered is a modification 
of the gear of case II. In case II the landing gear was considered to be 
connected to its supporting structure by a single springj in case III 
this single spring is replaced by a more complex structure described by 
some transfer function. This case is chosen mainly to demonstrate the 
application of the theory to complex problems for which the transfer- 
function concept may be of value. 


CASE I 

General Derivation 


In this section the basic equation of motion will be derived according 
to the summary theory for the special case of an inclined, rigid, swiveling 
landing gear (case I), which is illustrated in figure 5* This equation of 
motion could be obtained by making use of the previously derived equations 
of motion for the completely general case: however, it can more easily be 
derived here in a slightly different form for this particular problem. 
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The geometric quantities which enter the discussion of this particular 
landing gear are indicated in figure 5 • This gear has a swivel axis lying 
in the XY -plane and is inclined forward from the vertical Z-axis by a con- 
stant angle k . The perpendicular distance a between the ground-contact 
center point 0 and the swivel axis is called the trail. The swivel axis 
is assumed to move with constant velocity v along the X-axis without 
lateral motion from the XZ-plane. 

Rotation of the wheel structure about the inclined swivel axis by an 
amount i|r results in a component 0 of angular rotation about the verti- 
cal axis of magnitude 


0 = tjr cos k (ill) 

a component 7 of rotation about the X-axis (tilt) of magnitude 

7 = -t sin k ( 112) 

and a lateral deflection of magnitude 

r) 0 = -a\|r (113,) 

where all angles except k are considered small. 

The sum of all moments about the swivel axis must equal the inertia 
2 2 2 

reaction 1^0^^ = I^v D i|r, where 1^ is the moment of inertia of the 

wheel structure (including the wheel) about the swivel axis. The moments 
about the swivel axis are assumed to consist of the moments resulting 
from the previously discussed forces and moments that arise from tire 
distortion and ground loads plus the moments applied to the wheel by the 
supporting structure , which are assumed to consist of a restoring spring 
of moment pijr and a linear damper of moment g = gv Di|r> where p 
and g are constants. Thus, summation of the moments about the swivel 
axis gives the differential equation 


h Vo - , 'o) - V 


a - F. 


Vo + 1 1 - - v 


sin k + 


^(Djto - 9) cos k - rv 2 D(y 0 - t) 0 )cos k - pi|r - gv Df = I^D 2 ^ (114) 


where the first term is the total ground force due to tire lateral dis- 
tortion and tilt (see eqs. (20) and (27)) times its moment arm a; the 
second term is the vertical force times its moment -producing frac- 
tion sin k times its moment arm (see eq. (28)); the third term is the 
moment about the Z-axis due to tire twist (see eq. (24)) corrected by 
cos i|r for the component about the swivel axis; the remaining terms on 
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the left-hand side represent the gyroscopic torque due to lateral tire 
distortion (see eq. (32)), the spring restoring moment, and the linear 
damper moment. Now by making use of equations (ill) to ( 113 ), equa- 
tion (ll 4 ) can be written in the form 


Aj_ D i|r + Ar> Di|f + Aj\|r + DyQ + B 2 yQ = 0 


( 115 a) 


where 


A, = 


Ao = 


Az — 


B, = 


Bo = 


v 


aTY cos K + gv 

2x< r _l_ V si r* 2, 


a^K-^ + coshc + p + p K 
COS K + TV^COS K 


aK^ + CjF z sin k. 


(H 5 b) 


and 


p 

P K = aKy sin k - aF z sin k + ac-^F z sin k + c^F z sin k ( 115 c) 


The general relation between and yQ for this case is found by sub- 
stituting for tjq, 7, and 9 , according to equations (ill) to (113) in 
the general kinematic equation (19a). Thus 


V 2 - , . £Lh . . 

y Q + 2_ Z n D ^0 = " a + + l l^ cos K + r + sin K 

n=l 


or, abbreviating, 


and rearranging. 


ff = l + tan k 

r h 


^0Z^ cos k - a^\|r = 


n 


y 0 + z_ l n D y 0 = 2_ * n D y 0 

n=l n=0 


„n 


( 116 ) 


( 117 ) 


since Z^ = 1» Differentiating this result gives 


n+1 


{ol ± cos k - a^Dijr = l n D“'> 0 


n=0 
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which can also be written as 


OO 

cos K - a ) D t = XI i n_i D y 0 

n=l 


( 118 ) 


and, similarly, 

00 

(al ± cos k - a) D 2 \|r = XI l n -2 ( 1;L 9) 

n=2 

Substituting these relations into equation (115) and multiplying through 
by al - ^ cos k - a gives 


A 1 XI 

n=2 


n~2 


D n y r 




n 


dV 


+ 


cos k “ ajDy Q + ^V^l cos K “ aJy Q = ® 

Finally, after adding all terms of like order, substituting N = Z-^K^ 

(eq. (77 ) )) substituting for some of the A's, and using equation (ll6), 
there results the equation 

00 

H F n E n y 0 = 0 (120a) 

n=0 

where 


Fq = craN cos k + cos^k + p + p K + u^. 


F^ = a% + a'K^ cos k + p + P K ^q + gv + a2-jTv 2 cos 2 K 
F 2 = A x + A 2 Z^ + A jl 2 


(120b) 


F n = A l^n-2 + A 2 l n-1 + Vn 


(n > 2) 


and 


u K = c-^F z (ol^ cos K - a)sin k 


a 1 = a + (l - a)Z, cos k = a - sin k 


(120c) 
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Systematic Approximations 

Equations (120) provide the general differential equation of free 
motion for case I. The corresponding equations for the systematic approxi- 
mations A to D3 are obtained as follows. 

Approximation A .- The basic equation for approximation A is obtained 
by setting Z n equal to 0 for n > 3 in equations (120). The resulting 
equation can be stated in the following form: 


E o ^0 + E i D Vo + e 2 c5 yo + e 3 ® 2 yo + E 4 D yo + E 5yo = 0 


(121a) 


where 


E o = V 2 b 


E 1 = V 2 i 2 + (arv 2 cos k + gv)z^ 

Eg = + (arv 2 cos k + gv^Zg + ^a 2 K^ + cos 2 k + p + P K ^Z^ 

= I^v 2 + (aTv 2 cos k + gv)l 1 + ^a 2 K^ + cos 2 k + p + p K ^Zg 


E, 


2 p p 

E^ = a N + a'K a cos k + pZ-^ + p K Z 1 + gv + crl 1 TV cos k 


E^ = craN cos k + Kq, cos^k + p + p K + u r 
and 


> (121b) 


P K = (aK? - aE z + ac^F z + c^F z sin K^sin k 
u k = c^F z ^crZ.^ cos k - a^sin k 


(121c) 


Approximation B .- In order to obtain approximation B, set Z^ equal 
to 0 in equations (121b). This gives the differential equation 


E o + E i c3 yo 


+ Eg D y Q + E 5 Dy 0 + E^y Q = 0 


(122a) 


where 
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E 0 = 1 ^ v2l 2 

E 1 = I 2 . 1 + (arv 2 cos k + gv)z 2 

E 2 = I^v 2 + (arv 2 cos k + gv)l]_ + ^a 2 K^ + K a cos 2 k + p + P K ) i 2 
E^ = a% + a’K a cos k + pl-j_ + P K l l + + crZ-jTv^os 2 * 

E^ = aaN cos k + cos 2 k + p + p K + u^ 


(122b) 


(The corresponding equation with inclusion of hysteresis effects is 
listed in appendix B.) 

Approximation Cl .- The equation for approximation Cl is obtained by 
setting Zg equal to 0 in the equation for approximation B. The 

resulting differential equation is 


E q D 5 y Q + E 1 D 2 y Q + E 2 Dy Q + E^y Q = 0 


(123a) 


where 


E 0 = I^ 2 ^ 

E 1 = + ( aTv2cos K + gv)z-j_ 

E 2 = a 2 N + a'K a cos k + pZ-^ + p^Z 1 + gv + aZ^Tv^cos 2 * 
E^ = aaN cos k + cos 2 k + p + p k + u k 


(123b) 


Approximation C2 .- The equation for approximation C2 is obtained by 
setting | equal to 0 in the equation for approximation Cl. The 
resulting differential equation is 

Eq D^yg + Ej_ B 2 yQ + Eg DyQ + E^yQ = 0 ( 124 a) 


where 
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E 


0 


V 2 *1 




E, 


I^v 2 + (arv 2 cos k + gv) 


E, 


2 = a 2 N + al^ cos k + pZq + p^Zq + gv + Z-jTv^cos^k 


2 2 ., 


(124b) 


Ej = aN cos k + K a cos k + p + p K + u K q 


and 


u k1 = %*z{h cos K ~ a ) sin K (124c) 

(The corresponding equation, with hysteresis effects included, is listed 
in appendix B . ) 


Approximation Dl .- The differential equation for approximation D1 
is obtained by setting Z q equal to 0 in the equation for approxima- 
tion C2. The result is the differential equation 


E o ° 2 yo + E i D yo + E 2yo = 0 


( 125 a) 


where 


E r 


E-^ = a^N + al^ cos k + gv 




( 125 b) 


Eo = 


aN cos k + K cos^k + p + w, 

Ob n. 


and 


w^ = - aE z + c y& z sin *)sin k 


( 125 c) 


Approximation D2 .- The differential equation for approximation D2 
is obtained by dividing all terms in equations (124) by Z^ and then 

setting Zq equal to <» and using the relations N/Zq = and 

VZ X = €K^. The resulting equation is 
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where 


E, 


E, 


E, 


~r 

= v 

= arv 2 


cos k + gv 


(126b) 


= + aeK^ cos k + p + p K + tv 2 cos 2 k 


Ej = aK^ cos k + eK^ cosmic + c^F z sin k cos k 


Approximation D3 «- The differential equation for approximation D3 
is obtained by dividing the E's for approximation D1 by N and then 
setting N equal to <». The resulting differential equation is 


a Dy 0 + y Q cos k = 0 


(127) 


Previously Published Theories 

In the preceding section the differential equations are set down 
for case I according to the summary theory and the systematic approxi- 
mations thereto, complete with all pertinent details, including a number 
of second-order terms so as to enable the direct application of these 
equations to actual problems. In the present section the differential 
equations for the previously published theories are listed for the pur- 
pose of making clear the differences in the basic structure of the dif- 
ferential equations resulting from the application of these theories to 
case I. In order to avoid obscuring the more important differences 
between the equations of the various theories, all terms are omitted 
whose inclusion in any shimmy theory should be completely straight- 
forward (such as the spring re storing -moment coefficient p and the 
damping coefficient g) and also all inclination effects. Although the 
latter effects are not necessarily negligible, they do appear to be of 
second-order importance and their omission here should not alter the 
basic structure of these equations. With these omissions — that is, for 
g=p=K=0— the differential equations for case I according to the 
previously published theories are as follows in the terminology of the 
present paper: 

Von Schlippe -Dietrich and Rotta theories .- The basic equation of 
motion for this case according to the Von Schlippe-Dietrich theory, after 
neglect of the effects of tire width, corresponding to equation (120a) 
for the summary theory, is 
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where 


V 


V 


F ' 
n 


aN + 

,2 


aN + aK a + g^v 

V\-2 + V'n-l + ( a \ + K o,) l n* fc <*\(. h 


(n = 2, 


■ a ) 

• 0 


> (128b) 


F * = I , v^Z 0 + 
n V n-2 


g l vZ n-l + ( 


a Kv + K 1 - 

A ay n 


1 n-1 

ih Mli - a) 


(n - 3, 5, • • •) 


The symbol g-^ is a hysteresis damping coefficient defined by the equation 


g 


1 



/v 2 + d 2 2 co 2 


(128c) 


where d^ and d£ are hysteresis constants. Aside from the omission of 
some inclination effects and other terms, as discussed previously, the only 
differences between equations (128) and equations (120) for the summary 
theory lies in the inclusion of the hysteresis term involving g-^ and in 

the addition of the terms in brackets in equations (128b). Rotta's cor- 
responding equations, after neglect of tire width effects, would be the 
same as equations (128) except that Rotta omitted the hysteresis term. 

As was previously noted, the Von Schlippe -Dietrich theory and Rotta's 
theory differ only slightly in their respective considerations of the 
influences of tilt and tire width, neither of which effects are considered 
here . 


Bourcier de Carbon advanced theory .- The Bourcier de Carbon advanced 
theory leads to the fourth-order differential equation 


E o D Vo + E i D5 yo + e 2 c2 yo + e 3 D yo + E 4yo = 0 


(129a) 
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E 1 = V^l 

E 2 = I^v 2 + (a 2 K A + B^,)l 2 • 

E 5 = a% + aK^ 
e 4 - aN + ^ 


(129b) 


Bourcier de Carbon elementary theory .- The coefficients for the 
Bourcier de Carbon elementary theory are obtained by setting equal 

to 0 in equations (129) for the advanced theory. The resulting third- 
order differential equation is 

E 0 D 5y 0 + E 1 ^0 + E 2 Dy 0 + E 3 y 0 = ° (130a) 


where 


E o - V *1 

E 1 - v 2 

E2 = a% + al^ 


E3 = aN + K a 


(130b) 


Melzer .- Melzer's theory gives the third-order differential equation 


E 0 D?y 0 + E 1 c2y 0 + E 2 Dy 0 + E 3 y 0 = 0 


(131a) 


where 


E , 


E, = 


V 2i i 


(131b) 


E, 


2 = a% + 


E3 = aN + K a 
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Moreland advanced, theory .- Moreland's advanced theory leads to the 
fourth-order differential equation 


E 0 D Vo + E 1 E 5 y 0 + E 2 E 2 yo + E 3 Eyo + E 4yo = 0 


( 132 a) 


where 


E 0 * 

E i - Vi 1 ' 2 

E 2 = I^v 2 + a^TNv 
Ej = a 2 N + aTNv 


\ = aN 


(132b) 


Moreland elementary theory .- Moreland's elementary theory gives the 
differential equation 


a Dy 0 + yo = 0 


(133) 


Taylor . - Taylor’s theory leads to the third-order differential 
equation 

E 0 E 5 yo + E 1 E 2 y 0 + E 2 Dy 0 + E 3 y 0 = 0 


(134a) 


where 


E 0 = ^ 


E-l = 0 


E 0 = K^a' 


E, = La 
3 A 


(134b) 


Temple elementary theory .- The second-order differential equation for 
the Temple elementary theory is 


E q D y Q + E 1 Dy 0 + E 2 y 0 = 0 


(135a) 
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where 


E o = X V 2 

E]_ = a^N + aK a > 

Eg = aN + Kjj 


(135b) 


is 


Maier .- The second-order differential equation of Maier's theory 

-> 2 . 


E, 


'0 


d‘ y 0 + E i D y 0 + Vo = 0 


( 136 a) 


where 


E o = V 2 

E x = a% 
Eg = aN 


(136b) 


Kantrowitz .- Kantrowitz 
ential equation 

E o + E 1 


theory leads to the third-order differ- 
D 2 y Q + Eg Dy Q + E^y Q = 0 (l37a) 


where 

Eq = Iv L 
E^ = Iv 2 + aZgN + 
Eg — 0 
e 3 = aN + 


(137b) 


Wylie . - Wylie's theory leads to the third-order differential equation 
E q D 5 y Q + E-j_ D 2 y Q + Eg Dy Q + E^q = 0 (138a) 
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where 

E 0 = Iv 2 L 

E x = Iv 2 + aZ 2 N + Z 2 I^ 
E 2 = a% + aK^ 

E 3 = (aN + Ka) (l + 


Stability of Motion 

The basic equations of motion having been established according to 
the various theories for a rigid swiveling landing gear, attention will 
be directed next to the meaning of these different equations with respect 
to prediction of the shimmy behavior of the landing gear. However, before 
going into this subject in detail it may be useful to discuss briefly 
what sort of information is desired about the motion of a landing gear. 
Basically, the most important question is whether or not the motion is 
stable — that is, does the wheel tend to move in a straight line (with 
decaying shimmy oscillations or decaying aperiodic motion) or does the 
tire tend to move laterally out from its rectilinear course (with diver- 
gent shimmy oscillations or divergent aperiodic motion) . To answer this 
question of stability for linear systems, the analytic methods of Routh 
(ref. 26) or Hurwitz (ref. 2J) or graphical methods similar to those 
introduced by Nyquist (refs. 28 to 31) are available. A brief discussion 
of these methods is given in appendix C. Any of these methods will pro- 
vide for most cases a procedure for determining whether any particular 
combination of landing-gear parameters and rolling velocity is stable 
or unstable. 

In general, for complicated problems, rather than investigate the 
stability of a landing gear by these methods for all possible conditions 
it may be more convenient and sometimes more valuable to draw various 
types of stability diagrams describing the system in question. For 
example, for case I a typical experimental type of stability diagram is 
shown in figure J, which presents boundaries between the regions of sta- 
bility and instability as functions of trail and rolling velocity for a 
specific landing-gear model. Another useful type of stability diagram 
for some problems might be a plot of boundaries between stable and unsta- 
ble regions as functions of damping moment and rolling velocity. 

In determining these stability boundaries, use is made of the well- 
known fact that the motion of a linear system can change from a stable 
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to an unstable condition only where the motion is purely oscillatory — 
that is, where any variable, for example >(;, is of the form 

¥ = \ iVlX (139) 

or where the motion is purely uniform, of the form 

t = * m (i4o) 


Thus all possible stability boundaries can be obtained by directly sub- 
stituting expressions of the form of equations (139) and (lAo) into the 
basic differential equations. In connection with the question of what 
form of the differential equation to use, it is of some importance to note 
that the final form, where the equation is expressed in terms of one 
variable, is often not the most convenient form to use. For example, 
for case I the purely oscillatory boundaries are most advantageously 
obtained by using equations (115) and (117) with the substitutions 


* = V 


iVqX 


A 


'0 


i (viX+aJ 

y n e x/ 


(141) 


The advantage of this particular choice arises from the fact that it leads 
to two algebraic equations, one of which does not include the damping 
parameter g. This isolation of the parameter g usually eases slightly 
the mathematical labor of solving for the purely oscillatory boundaries. 

The equations governing the stability boundaries for case I for the 
general theory and for the systematic approximations are listed in 
appendix D. 


Comparison and Evaluation of the Summary Theory and Its 
Systematic Approximations 

The dual object of the present section is (l) to assess further the 
value of the summary theory by comparisons between the predictions of 
this theory and the available experimental data for case I conditions 
and (2) to determine, by comparison of the relative predictions of the 
summary theory and its systematic approximations, what is the simplest 
satisfactory systematic approximation to the summary theory. Discussion 
of the previously published theories, as applied to case I conditions, 
is contained in a subsequent section. 
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For convenience the following discussion Is divided into separate 
considerations of stability-boundary conditions and unstable shimmy 
conditions. 

Stability-boundary condit ions . - The present section discusses theo- 
retical and experimental stability-boundary conditions insofar as they 
are influenced by the tire parameters Z n (where n = 1, 2, . . .), 4 , 

K / e A , N, and t and by hysteresis effects. In the major part of 

this discussion the stability boundaries to be considered will be of the 
type . obtained by plotting curves of trail against rolling velocity for 
those trail conditions that separate regions of stability and instability. 
The general shapes of these stability boundaries for case I, according 
to the summary theory and the systematic approximation theories A to D3, 
are sketched in figure 8 for the special condition of no damping or gyro- 
scopic moments (g = r = 0) . It is seen that the summary theory and 
approximations A to C2 each predict that at high speeds the motion is 
stable for large trails and unstable for small trails; the horizontal 
boundary line is the same for each case, and is generally located at a 
trail roughly equal to the tire radius. (This boundary is theoretically 
completely independent of the spring restoring moment p D^-ijr and is 

relatively independent of swivel-axis inclination k.) Approximations Dl, 
D2, and D3 fail to predict this boundary. Also, these three approxima- 
tions, together with approximations Cl and C2, fail to predict any effect 
of rolling velocity on the low-speed stability boundaries, whereas, 
according to approximation B, for sufficiently small speeds the motion 
becomes stable for all small trails and, according to the higher theories, 
for most of the small-trail region. Also, at low speeds and large (usually 
impractical) trails the higher theories (B and above) indicate that the 
motion becomes unstable at sufficiently small speeds. The effects of the 
omitted damper and gyroscopic -moment terms would be to reduce the size of 
the regions of instability. 

(a) Effect of higher Z n terms: As a first test of the summary 

theory and its systematic approximations there are available the experi- 
mental data of Von Schlippe and Dietrich (refs. 3, and 5) which were 
obtained with a small model landing gear equipped with a pneumatic tire 
of 26 cm (10 in.) diameter. This model landing gear was tested at rela- 
tively low speed conditions where the higher Z n terms (z 2, Zj, . . 

are of some importance; consequently, these data provide an opportunity 
for testing the relative and absolute validity of the summary theory and 
the higher approximations A to C2 (which differ essentially only by their 
inclusion or omission of the higher Z n terms). 

The basic landing gear and tire constants for the Von Schlippe- 
Dietrich model, which was tested only in the untilted condition (k =0), 
as taken from references 3 and are as follows: 
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K = P = g = 0 

1-^ « 0.53 + 0.0025a 2 cm-kg-sec^ 

L = 10 cm 
N = 64 0 kg/radian 
= 3 j> 040 cm- kg/radian 
= 45 kg/cm 

The quantities l-,, h, and the higher l n 's were calculated from the 
previously discussed relations Ip = n/K^, h = Zp - L, and 
l n = (nL + h)h n "^/nl (See eqs. ( 77 ) and (19a).) 

The experimental data obtained by Von Schlippe and Dietrich for 
this model are shown in figures 9 and 10, together with the corresponding 
predictions of the summary theory and the systematic approximations A 
to C2. (Also shown are the predictions of the theory of Von Schlippe 
and Dietrich which are discussed in a later section.) Figure 9 presents 
stability-boundary plots of trail against velocity, and figure 10 pre- 
sents the frequency at these stability boundaries as a function of veloc- 
ity. No theoretical curves are shown on these figures for approxima- 
tions Dl, D2, and D3 since these approximations are too crude to give 
any detailed information for this problem; they predict either completely 
stable or completely unstable motion for all positive trails (see fig. 8). 
The equations used to calculate the theoretical curves in these two fig- 
ures are given in appendix D. In these calculations the gyroscopic torque 
term involving t has been neglected since t is unknown for these data. 
A rough value of t could perhaps be estimated, but such a dubious esti- 
mate did not appear necessary because the term involving t, according 
to any reasonable estimate of t, would be of no importance in the veloc- 
ity range of these experimental data. 

In comparing the theoretical curves in figures 9 and 10 it is 
observed that approximation A gives a boundary very close to that of the 
s umma ry theory. Approximation B does not give as close agreement but it 
is still fairly good and, more important, for most of the trail range 
the difference between approximation B and the summary theory is small 
as compared with the difference between the summary theory and the experi- 
mental data. As was previously noted, approximations Cl and C2 (which are 
identical for the present condition of k = 0) predict a trail- velocity 
stability boundary which is independent of velocity so that this approxi- 
mation is an inadequate representation of the summary theory at low veloc- 
ities. However, at high speeds approximations Cl and C2 give the same 
stability boundary and frequency as the higher approximations. 
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As a further aid in comparing the different systematic approximations 
with the summary theory, figure 11 presents a plot of the linear damping 
coefficient g required to stabilize the motion of the Von Schlippe- 
Dietrich model at a medium trail of 7 cm as calculated according to the 
summary theory and the various systematic approximations (the equations 
used are presented in appendix D) . This figure confirms the conclusions 
drawn from figures 9 and 10 that approximation A is a very good repre- 
sentation of the summary theory and that approximation B is also a good 
representation of the summary theory. However, more importantly, fig- 
ure 11 demonstrates that approximations Cl and C2 also give a fairly good 
representation of the summary theory with respect to prediction of the 
maximum amount of damping (i.e., the maximum value of g) required for 
stabilizing the motion. Approximations Dl, D2, and D3 are seen to give 
inadequate representations of the summary theory. 

The preceding conclusions are, of course, only proven to be valid 
for the specific conditions of the Von Schlippe-Dietrich model tests. 
However, they are believed to be valid for most practical rolling 
conditions . 

In considering the correlation between theory and experiment for the 
Von Schlippe-Dietrich test conditions, it is noted that the experimental 
stability boundary in figure 9 is of the same general shape as that given 
by the summary theory and approximations A and B but that it lies to the 
right of the theoretical curves and thus indicates that the experimental 
system is more stable than the theoretical system. Similarly, the experi- 
mental frequency- velocity curve in figure 10 falls below the theoretical 
curves. These discrepancies are perhaps a result of the neglect of hys- 
teresis damping in the calculation of the theoretical curves of figures 9 
and 10 and are discussed more fully in the next section. 

(b) Effect of hysteresis: In order to investigate whether the dis- 

crepancies between theory and experiment shown in figures 9 and 10 might 
be explained by a consideration of hysteresis effects, some results of 
calculations involving hysteresis effects are shown in figures 12 and 13 
for the same test conditions as in figures 9 and 10. Figures 12 and 13 
present theoretical calculations of trail-velocity and frequency-velocity 
stability boundaries, respectively, both with and without consideration 
of hysteresis effects, together with Von Schlippe and Dietrich's experi- 
mental data. The curves for the Von Schlippe-Dietrich theory were cal- 
culated by using Von Schlippe and Dietrich's theory with their hysteresis 
equations (see eqs. 128), whereas the curves for approximation B were 
calculated by using the equations for hysteresis effects derived in the 
present paper (see appendix B). The values of and indicated 

in these figures, which were used in these calculations, were estimated 
from static hyst ere sis -loop data, partly unpublished and partly published 
in reference Zk. The values of dj and d 2 used were taken from refer- 
ence 5. It is observed from these figures that there is a considerable 
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difference between the calculations according to either theory, depending 
on whether hysteresis effects are included or omitted. Moreover, it is 
seen that either theory provides fairly good agreement with the experimen- 
tal data when the hysteresis effects are taken into account. These facts 
suggest that hysteresis effects may be an important factor which must be 
taken into account in order to obtain good agreement between theory and 
experiment. On the other hand, neither the Von Schlippe -Dietrich nor the 
present paper's consideration of the hysteresis effects seems to rest on 
a completely sound foundation. Thus it appears safe to conclude only that 
since two different approaches to the hysteresis problem indicate that 
hysteresis effects are important, a more rigorous analysis of the hysteresis 
problem would be worthwhile. 

(c) Effect of Zp: The next test of the summary theory makes use of 

the experimental data of Melzer (ref. 10), who performed a series of tests 
with an untilted (k = 0) solid rubber tire 7 cm. (3 in.) in diameter at 
sufficiently- high speeds so that his data would be expected to fall con- 
siderably to the right of the curved- line low-speed stability-boundary 
curve in the first two sketches of figure 8. For this velocity range the 
predictions of the summary theory and approximations A to C2 are identical; 
they all predict that the undamped (g = t = T^ = T a = 0) motion is stable 
for all values of trail a greater than Zp = e/K^; that is, the critical 
value of trail a c at which the motion changes from an unstable state to 
a stable one is given by the relation a c = Zp. (Even if gyroscopic moments 

are taken into account, this equation is only slightly modified throughout 
a relatively large range of rolling velocity; however, for very large veloc- 
ities this relation breaks down as a result of the gyroscopic moment and 
the motion becomes stable for all positive trails. This phenomenon will 
be discussed in a subsequent section.) The prediction a c = Z-^ is well 

confirmed by Melzer' s tests, as is illustrated in the following table of 
data taken from reference 10, which lists the experimental values of Zp 

together with the trail required for stability for several conditions of 
vertical loading: 


F z , kg 

1 

2 

2.8 

3.6 

h = N / K x> cm 

2.58 

2.85 

3.22 

3.40 

a c , cm 

2.5 

3.0 

3.1 

3.4 


These particular data were taken for the case where no spring restoring 
force acted on the model landing gear (p = 0). Similar good agreement 
was obtained for the case where a strong linear spring restoring moment 
was present (p > KqJ . For this case, according to the summary theory 
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and approximations A to C2, the stability boundary for the positive trail 
condition is the same as for the case of no spring restoring force. This 
prediction is well confirmed by Melzer's tests, as is illustrated in the 
following table of data also taken from reference 10. 


F z , kg 

1 

2 

2.8 

3.6 

Z x = N/K^, cm 

1.96 

2.60 

2.90 

2.91 

a c , cm 

2.0 

2 - 2.25 

«2.7 

«2.8 


(The difference between the values of Ip in these two tables is merely 
a consequence of changes in tire characteristics between the corresponding 
tests. ) 

In order to assess the significance of the preceding comparisons of 
theory and experiment it should be noted that the theoretical relation 
a c = l 1} calculated for g = r = T^ = T a = 0, is independent of spring 

restoring moment p and tire torsional stiffness I 1 ^ and, in the veloc- 
ity range discussed, can be shown to be not strongly influenced by gyro- 
scopic or hysteresis moments (t, T-^, and T a ) or by the higher Z n 's 

(Z 2 , Z^., . . .)• Consequently, this comparison tells practically nothing 

about the correctness of the manner in which these important quantities 
have been inserted into the summary theory. On the other hand, the theo- 
retical relation a c = Zp depends almost entirely on the correctness of 
the kinematic equation of the summary theory which, for Melzer's test 
conditions, reduces to the kinematic equation of approximation C2 
(eq. (85)): 


Aq — — Z x a 

Thus, the results of the preceding comparison indicate that there exists 
a range of rolling speeds in which the kinematic equation of the summary 
theory, as well as of approximations A to C2, is reasonably correct 
(except possibly for the terms involving £, which are as yet not evalu- 
ated and are not very important). 

It can be said with safety that the range of velocity for which the 
theory gives good agreement with Melzer's model data corresponds to full- 
scale conditions somewhere inside the practical rolling speed range and 
possibly covering much of the practical range. However, the preceding 
comparison definitely does not prove anything about the adequacy of the 
s umm ary theory for small velocities or for the highest velocities which 
may be encountered in practice. 
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Further confirmation of the preceding conclusions is provided by 
the experimental data of Schrode (ref. 17) who performed tests, similar 
to those of Melzer, for realistic pneumatic tires as large as 39 cm (15 in. 
in diameter, as compared with the small 7 cm (3 in.) solid rubber tire 
tested by Melzer. Schrode obtained trail-velocity stability-boundary 
plots of the type illustrated in figure 7* These stability -boundary plots 
indicate the same result as Melzer' s data, namely, that there exists a 
range of velocity in which the motion is stable above a certain critical 
trail a c and unstable below it. (This velocity range for the data in 
figure 7 is approximately 60 to l60 km/hr.) It is not possible to check 
quantitatively the theoretical stability-boundary equation a c = Z-^ for 

Schrode 's data because Schrode provides no information suitable for accu- 
rately evaluating Ip. However, some qualitative confirmation may be 
found, since the quantity Zp always appears to be of the order of magni- 
tude of the tire radius r and for Schrode' s data a c is found to be of 
this same order of magnitude (e.g., see fig. 7)* Thus Schrode' s experi- 
mental data appear to confirm the previously drawn conclusion that there 
exists a velocity range in which the kinematic equations of the summary 
theory and approximations A to C2 are valid. 

Dietz and Harling have presented some similar stability-boundary 
curves in reference 16 which also confirm the foregoing conclusions. 

(d) Effect of *■ : Some insight into the effect of the tilt param- 

eter | can be obtained by an examination of the effects of swivel-axis 
inclination k. on the stability boundaries according to the predictions 
of approximation Cl for the condition where damping, spring-restoring, 
and gyroscopic moments are neglected (g=p=T=0) in order to isolate 
the effects of inclination. (These assumptions appear to be justified 
for the experimental conditions discussed in this section.) Under these 
assumptions one theoretical stability boundary is given by the equation 

a c = Zp cos k + - ^ sin tc (l^2) 

Experimental data suitable for testing this relation are available in 
reference l6 for an inclination range of -20° < k < 20° for one constant - 
velocity condition. These experimental data, some of which had to be 
slightly extrapolated from the data in reference 16 , are shown in fig- 
ure it together with the predictions of equation (lt2) for values of | 
equal to 0 and 1. Although reference l6 does not supply the values of L, 
h, and Z^ needed for calculations, the assumed values indicated on the 

figure are probably accurate enough to justify the following more or less 
qualitative observations. (The value of Zp was chosen so as to make the 
calculated and experimental values agree for the case k = 0.) It is noted 
that the experimental variations and the theoretical variations for | = 0 
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are in fairly good agreement and also that these two variations are more 
or less symmetrical with respect to positive and negative values of k» 

On the other hand, the theoretical curves for | > 0, such as the indicated 
curve for | = 1, will all he unsymmetrical. Thus, it appears that | is 
probably close to zero. In this connection it might be noted that 
Greidanus' theory, which is the only known theory that uses a term cor- 
responding to I, implies a value of i > 1 (compare eqs. (100) and (101)). 

(e) Effects of Ky and c-^i In order partly to assess the importance 
of the tire parameters and c^, the special case of approximation C2 

where r = p = 0 and l ^ cos k - a = 0 will be considered. While this 

particular case is of no practical importance in itself, its examination 
permits some insight into the effects of the tire parameters K ^ and c^. 

For this case the damping required to stabilize the motion is given by the 
inequality 


a(aN + cos k 'j + P K l] 


gvi ± + i^gv 5 + g 2 y 2 ^ 1 > 0 


where is given by equation (115c). It is evident from this relation 

that p„ is a stabilizing term if positive and destabilizing if negative. 

ft 

Also, according to reference 2, may be as large as F z and, according 

to reference 2b, c-^ O.75. Thus, according to equation (115c), if the 

small sin 2 K term is neglected, p K may be positive for positive k, 
whereas if the K^, 0^, and c terms are neglected (as has been done 

in all previous investigations except ref. 2) then P K is always negative 
for positive k. Therefore, if situations should arise wherein p K is 
important, it is not necessarily safe to neglect the terms involving c^ 
and K^, that are used in the determination of p^ . (See eq. (115c).) 

(f)- Effect of cornering power N: As a rough check on the variation 

of the tire cornering power N under dynamic conditions, there are avail- 
able experimental frequency data obtained by Melzer in connection with his 
previously mentioned tests on an uninclined (k = 0) model landing gear 
equipped with a solid rubber tire of 7 cm (3 in.) diameter (ref. 10). For 
the higher velocity conditions of Melzer* s tests, the predictions of the 
summary theory and approximations A to D1 lead to the frequency equation 


■f - 1 , IbN + Kq + p 

2*1 1 it 


(1^3) 


for an uninclined and undamped landing gear, that is, for k = t = g = 0. 
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(inclusion of the effect of finite t in this equation would not signif- 
icantly alter this equation for the test conditions to be discussed herein*) 
Some of Melzer's experimental data are compared with the predictions of 
this equation in the following table for the condition p = 0. The experi- 
mental data shown represent Melzer's data for the highest velocity condi- 
tions tested. 



2.8 

3.6 

a /h 

0.47 

O .78 

0.44 

0.73 

0.88 

^calculated' C P S 

3-8 

4.5 

4.0 

4.8 

5.1 

^experimental' C P S 

3-3 

3-5 

2.7 

4.1 

4.7 


The theoretical and experimental values are seen to be in fair agreement. 
However, the experimental values are somewhat smaller than the corresponding 
theoretical values. This discrepancy is believed to be largely due to the 
fact that these experimental tests were not conducted at sufficiently small 
values of shimmy amplitude for the assumptions of a linearized theory to 
be valid. Specifically, all of Melzer’s frequency data were obtained for 
maximum swivel angles of 5° or larger. (The data shown in the preceding 
table correspond to the condition 8 m = 5°*) Moreover, Melzer's data 

indicate that there is a fairly definite decrease in shimmy frequency with 
increasing maximum swivel angle. A sample plot of Melzer's data illus- 
trating this effect is given in figure 15. Also shown is the theoretical 
calculation, which is valid only for 8 m = 0°. If allowance is made for 
a certain amount of experimental error, extrapolation of the experimental 
data to 8 m = 0° could be considered as confirmation of the theory. It 

should be noted, however, that much of the rest of Melzer's data, while 
not necessarily disputing this conclusion, do not so clearly support it. 
Also, plots of the type of figure 15 are of limited significance since 
each test point shown corresponds to a different rolling velocity. In 
view of these considerations, the only reasonable conclusion that can be 
reached appears to.be that Melzer's data roughly confirms the theoretical 
frequency and does not conclusively dispute its quantitative accuracy. 

Melzer has also conducted frequency tests on the same model with an 
additional strong restoring spring (spring stiffness several times the 
tire torsional stiffness). A comparison of theoretical and experimental 
frequencies for this test is shown in the following table: 
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F z , kg 

2.0 

2.8 

3.6 

a / l l 

0.77 

O .69 

0.86 

0.69 

0.86 

^calculated 3 9P S 

5-2 


5-7 

5.5 

5.8 

experimental* 9P S 

^•9 

5^5 

5*9 

5*8 

5*85 


The much better agreement obtained for this case is explained by the pre- 
dominant influence of the spring restoring moment, since for large values 
of p the model system approaches the condition of a simple torsional 
oscillator with moment of inertia 1^ and spring constant p, for which 

condition the well-known frequency equation is 2rtf = 

In order to assess the significance of the preceding comparisons, 
the quant it ies involved in the theoretical equation (lV5) — a, N, 

p, and 1^ — will be considered. The quantities a, p, and 1^ are 

easily measured constants and, for most of Melzer's data, is much 

smaller than aN; therefore, the preceding fair agreement between theory 
and experiment indicates that the quantity N (the tire cornering power) 
which was considered to be a constant in the preceding calculations, 
actually does not vary extremely with rolling velocity and shimmy fre- 
quency — at least not for Melzer's test conditions. 

(g) Effect of gyroscopic torque: The next question to be considered 

is the influence of the gyroscopic torque resulting from lateral distor- 
tion of the tire. All pertinent experimental data obtained at very high 
speeds (e.g., ref. 10; see also fig. 7) demonstrate that at sufficiently 
high speeds the previously discussed conclusion that high-speed motion is 
unstable for trails less than is no longer valid. Instead, the 

experimental data show that at these very high speeds instability at any 
given positive trail ceases above a certain critical velocity. The exist 
ence of this critical velocity will now be shown to result, at least in 
part, from the gyroscopic action which was previously included only in 
Kant rowitz 1 theory (ref. 8) but was not specifically mentioned there. 

The simplest systematic approximation that adequately provides for this 
effect is approximation C2. In order to isolate the gyroscopic effect, 
the special condition of no tilt (k = 0) and no spring-restoring force 
(p = 0) or damper (g = 0) will be considered. For this condition the 
equation for the stability boundary of approximation C2 (or Cl) reads 

(Vc 2 + aTV c 2z i )( a% + aK a + ¥v) = Vc^lt 8 ® + K a) 
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where the underlined terms are the gyroscopic terms. For the computation 
of the critical velocity v c this equation may he simplified still further 

if it is realized that the quantity arl^ is small in comparison with the 
moment of inertia 1^ about the swivel axis; hence , for an approximate 
calculation the term aTv 2 !^ can omitted. Then solution of equa- 
tion (1^4) for the critical velocity v c above which the system is stable 
yields the expression 

v 

which is observed to give an infinite critical velocity for zero gyroscopic 
action (t = 0). 

The only available experimental data containing enough information 
on the tire constants that are necessary for checking the validity of 
equation (l45) are Melzer's data (ref. 10) and even these data do not 
provide the required gyroscopic moment; therefore, it can only be crudely 
estimated as follows: The mass of the tire will be of the order of magni- 

tude W 1 j^2jt (r - r^Jjitr^ 2 , where r is the overall tire radius, r^ the 
radius of the cross section of the tire torus, and w^ the average tire 

density. The moment of inertia will be the mass times the radius of 
gyration r g squared; thus, with = 1/2 according to Kantrowitz, 

t (eq. (33)) becomes 

T = - r Q r g 2 

r(r + r^) 

For the usual tire r^ 0.3r, r^ is slightly smaller than r (say 
r ^ ~ 0 . 9r ) , and r g is probably around 0.8r. Then, to a crude approxi- 
mation, t « 0.21w-jjr5. For Melzer's solid rubber tire r = 3°5 cm and 
w 1 is probably about 10 kg-sec 2 /cm^ (specific gravity of 1); thus 

t «= 10"5 kg-sec 2 /cm. Critical velocities calculated from equation (1^5) 
with this value of t are compared in figure 16 with some of Melzer's 
experimental data for one test condition at various values of a /l]_* The 

calculated and experimental values of critical velocity are seen to be of 
the same order of magnitude. Since neglect of the gyroscopic moment 
gives, theoretically, an infinite critical velocity, this agreement indi- 
cates that the gyroscopic moment is an important factor in producing sta- 
bility at high velocities. Also, the theoretical calculation is conserva- 
tive, that is, the unstable region is overestimated. The quantitative 
agreement between theory and experiment is fair but far from excellent. 
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One probable reason for some of the disagreement is the relatively crude 
procedure used for estimating the parameter t. Another possible explana- 
tion may lie in hysteresis effects, as follows; 

If the differential equation for approximation C2 (see eq. (124a)) 
is modified to take into account hysteresis forces and moments in the 
manner suggested in this paper, a modified differential equation results 
(see eq. (B2)), which has the stability boundary equation 

(Vc 2 + arv 2 ^ + a 2 NT A v c + aK^T^ ^aSj + aK^ + 

Z lT v c 2 + aNT A v c + K^Tv^ = + K^,) 

for the same conditions as the corresponding equation (l44), 

K = p = g = 0. After neglecting the underlined terms, which 
small at large velocities, this equation can be expressed in 
and more easily interpreted form 

2 + + y, _ (l ! - a)(aH + ) 

C Z 1 T C Z 1 T 

This equation indicates that the effect of finite hysteresis 
T-^ 4 - 0; T a ^ oj is to reduce the critical velocity below what it would 

be for no hysteresis effect (t^ = T a = 0). This result is also indicated 

in figure 16, where calculated and experimental curves are shown for the 
previously discussed high-speed conditions of Melzer's model tests. (The 
values of T^ and T a needed for these calculations were obtained from 

equations (60), (62), and (l43), by using the previously mentioned esti- 
mated values of = 0.1 and = 0.2, based on static hysteresis 

loops . ) In figure l6 the experimental data lie between the theoretical 
curves for "hysteresis considered" and "hysteresis not considered." The 
theoretical calculation that includes hysteresis is extremely unconserva- 
tive. Two conclusions can be drawn from these observations. First, if 
the actual hysteresis effect at high speeds is only a fraction of the cal- 
culated effect, this fact might explain the difference between the experi- 
mental curve and the theoretical curve that does not include hysteresis. 
Second, as was previously noted, it is evident that the treatment of 
hysteresis effects in the present paper is inadequate and unconservative 
at high velocities. 

In concluding this discussion of gyroscopic torque, it should be 
noted that for the case of a rigid landing gear the critical design con- 
dition (velocity at which sh immy is most intense) occurs at low rolling 
speeds where the gyroscopic moment is insignificant. Thus, the inclusion 


(146) 

namely, 

are relatively 
the simpler 


(147) 

(that is. 
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of this gyroscopic moment in the theory is of somewhat academic interest 
(at least for case I) and it probably could be safely omitted in practical 
design calculations. 

Unstable shimmy conditions .- As a further overall check of the sum- 
mary theory and its systematic approximations, the experimental data of 
Kantrowitz (ref. 8) for unsteady shimmy conditions are available. 

The significant features of unsteady shimmy motion are the divergence 
and frequency of the oscillation, which are simply the real and imaginary 
parts of the roots of the characteristic algebraic equation corresponding 
to the differential equation in question. Kantrowitz has made measure- 
ments of these quantities for a model tire of 4-inch diameter at inclina- 
tion angles k of 5° and 20° with corresponding trails of about 0.08r 
and 0.31r, respectively (ref. 8). His experimental results for k. = 5° 
are presented in figure 17, together with corresponding theoretical cal- 
culations made according to approximation B, which is the simplest system- 
atic approximation to the summary theory which describes, at least quali- 
tatively, the sh immy phenomena throughout the complete range of rolling 
velocity. The theoretical and experimental frequencies are seen to be in 
fairly good agreement. The theoretical and experimental divergences are 
in fair qualitative agreement, but the experimental variation is sometimes 
considerably below the corresponding theoretical one. This quantitative 
disagreement may be due to several factors. First, hysteresis effects 
are neglected in the theoretical calculations. Although use of the hys- 
teresis force and moment equations derived in this paper would not com- 
pletely explain the disagreement, it is believed that these hysteresis 
equations are not accurate enough, particularly at small trails (a = 0.08r 
for the data in fig. 17), to justify the conclusion that the disagreement 
cannot be explained by hysteresis effects. A second partial explanation 
of the disagreement arises from the fact that the theoretical calculations 
may be based on insufficiently accurate values of the necessary tire param- 
eters, since Kantrowitz did not make direct measurements of all of the most 
fundamental tire parameters, such as h, a, N, and 1^. Specifically, he 

measured only the quantity L, a quantity approximately equal to 

aN cos k + cos^ k for 2 values of k, and the path frequency Vp and 

trail a for kinematic sh immy (shimmy with velocity approaching zero). 

The basic tire parameters used for calculating the theoretical curves in 
figure 17 were approximately deduced from these quantities as follows: 

The quantity h was obtained from equation (Dl) after setting v equal 
to 0 and substituting Kantrowitz' experimental values of L, Vp, and a 

for kinematic shimmy. This procedure for determining the quantity h is 
not necessarily accurate since equation (Dl) neglects hysteresis effects. 
The tire deflection needed for calculating the trail was estimated from 
figure 8 of reference 2. The trail was computed from the tire radius, 
the tire deflection, and the inclination. With the aid of this estimated 
value of trail the tire parameters N and 'K^ can be obtained from 
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Kantrowitz' approximate values for aN cos k + cos k. In addition, 
most tire tilt effects were neglected) specifically, £ and p K were 

taken equal to zero. While the foregoing procedure will probahly give 
roughly correct values for most of the fundamental tire constants, it is 
believed that the limitations of this procedure and the neglect of hys- 
teresis effects are sufficient to prohibit the making of any strong point 
out of the discrepancies between theory and experiment in figure 17* Thus, 
to summarize, it appears that Kantrowitz' data furnishes only a rough over- 
all confirmation of the summary theory. Although quantitative agreement 
is poorer than for most of the previously discussed experimental data, 
this poorer agreement is not necessarily significant. 

This completes the discussion of case I with respect to the summary 
theory and its systematic approximations. Next, attention will be directed 
to a discussion of case I with respect to the predictions of some of the 
previously published theories. 


Predictions of Some of the Previously Published Theories 

Some interesting features of the previously published theories in 
relation to case I are discussed in the following paragraphs. Comments 
on the influence of swivel-axis inclination will not be repeated 
here. 


The theory of Von Schlippe and Dietrich gives predictions which are 
substantially the same as the predictions of the summary theory, as can 
be seen by a comparison of the predictions of these two theories in fig- 
ures 9 and 10 for Von Schlippe and Dietrich's model test conditions. 
Hysteresis effects were neglected in computing both sets of theoretical 
curves. In comparing these two theories it should be noted that the only 
difference in the two sets of theoretical curves arises from a slight dif- 
ference in the expressions used for the elastic forces and moments of the 
tire (see section entitled "Forces and Moments on the Wheel" and the com- 
ments after equation (128)). The Von Schlippe-Dietrich theory also pro- 
vides for some tire width effects, but these effects are believed to be 
relatively small for the present test conditions and were not taken into 
account in computing the theoretical curves In figures 9 and 10. From 
these figures it is seen that the differences between the stability bound- 
aries and frequency curves for the Von Schlippe-Dietrich theory and the 
s ummar y theory are usually small in comparison with the differences between 
the theoretical curves and the experimental data. Thus, it seems reason- 
able to conclude that there is no significant difference between the main 
features of the summary theory and the Von Schlippe-Dietrich theory, 

Bourcier de Carbon's advanced theory provides essentially the same 
predictions as approximation B and will thus probably give a reasonable 
prediction of sh imm y behavior for the complete velocity range. Similarly, 
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Bourcier de Carbon's elementary theory, corresponding to approximation C2, 
will probably give reasonable predictions for the high velocity range. 


Melzer's theory correctly predicts the existence of the large-trail 
stability boundary given by the equation & c = but it also predicts 

the existence of stable motion in the small negative trail region between 
zero trail and a trail, equal to -e = -K^N. The latter prediction is in 

disagreement with the experimental data of Von Schlippe and Dietrich 
(ref. 5) who conducted some tests in this trail range and found the motion 
there to be unstable. 


The stability boundary according to Moreland's advanced theory for 
the case of no damping or spring-restoring forces (see eq. (l32a)) is 
given by the equation 


where 



NZ-jT 2 



(JA8) 


This equation is plotted in figure 18 for zero time constant (for which 
case Moreland' s theory reduces to the subcase of approximation C2 for 
which e = Kq^ = 0) and for several finite values of the time-constant 
parameter T2« It is seen that if Tg is large, a large-trail stability 

boundary no longer exists at the trail a c = Since the actual exist- 

ence of this large-trail stability boundary has already been demonstrated, 
it appears likely that Tg cannot be very large. On the other hand, if 

Tg is small the introduction of the time-constant term produces an almost 

linear decrease of critical trail with increasing velocity until a certain 
limiting velocity (equal to lq/T) is reached^ above this velocity all 
motion is stable. Thus, the influence of the time-lag constant term is 
somewhat like that of the previously discussed gyroscopic moment due to 
tire distortion or the hysteresis force and moment, which may also produce 
stability at high velocities. However, in regard to the general shape of 
the critical trail-velocity curve, the variations predicted by considera- 
tion of the gyroscopic or hysteresis effects (see solid and dashed lines 
in fig. l6) appear more like those of the experimental data (figs. 7 or 16) 
than does the nearly linear variation predicted from Moreland's time-lag 
term for small T2» 


Moreland's elementary theory, Temple's elementary theory, and Maier's 
and Taylor's theories are too crude to give any details for case I. 
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Kantrowitz’ theory incorrectly predicts instability for all positive 
trails in the absence of damping or gyroscopic moments. (This prediction 
is a consequence of the fact that the third coefficient in Kantrowitz' 
differential equation ( 137a) is zero.) 


Wylie's theory (see eqs. (138)) correctly predicts the existence of 
stability at large trails; however, the predicted value of critical trail 
is given by the equation 


a c (a c + e)NZ 2 = Iv 2 L (l49) 

for k = 0. This relation implies that the critical trail is a continu- 
ously increasing function of velocity, whereas the previously discussed 
experimental data clearly indicate that the critical trail rapidly reaches 
the maximum value Z-j_. Also it is noted that the last of Wylie's equa- 
tions (138b) contains the factor 1 + a/L, which does not appear in any 
of the other theories. This term appears as a consequence of the earlier 
mentioned fact that Wylie's theory does not correctly predict the influ- 
ence of trail for the yawed rolling condition. It might be noted that if 
Wylie's theory were modified to remove this difficulty, as suggested in 
equation (llO), then this factor 1 + a/L would be replaced by 1; thus 
the modified Wylie theory would be more consistent with the other equations 
Also, for high velocities the modified Wylie theory would more correctly 
predict a definite critical trail according to the relation a c = L. 


Practical Application 

Before concluding the discussion of case I, it is perhaps pertinent 
to comment on whether the preceding theoretical predictions for this 
idealized case can be applied to practical landing-gear problems. Moreland 
has demonstrated that neglect of the torsional and lateral elasticity of 
the landing-gear strut can sometimes lead to false predictions of stability 
in an actually unstable system (ref. 12). For example, for too large tor- 
sional damping the torsional damper unit effectively locks the swiveling 
structure against torsional movement with respect to its point of attach- 
ment to the landing gear strut, so that torsional movement of the swiveling 
structure is possible only as a consequence of the always finite torsional 
stiffness of the strut. In other words, for too large damping the damper 
unit and strut combination behaves like a pure torsion spring and is thus 
ineffective for damping purposes. In regard to lateral strut stiffness 
Moreland has pointed out that for a rigid tire, corresponding to approxi- 
mation D3, if the strut is prohibited from lateral motion, as is assumed 
for case I, the motion is stable; on the other hand, if the strut is 
assumed to be of stiffness approaching infinity, the system may be unsta- 
ble, While this particular criticism of the practical value of case I 
for infinite stiffness applies only to approximation D3, still for finite 
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strut stiffness of the same order of magnitude as the tire lateral stiff- 
ness the theoretical equations of the other approximate theories may also 
be unconservative. 


CASE II 


This section is concerned with the idealized landing gear shown in 
figure 6. This landing gear consists of a wheel free to swivel about an 
uninclined always-vertical swivel axis that is attached by a horizontal 
linear spring, of spring constant k, to the supporting structure. This 
configuration, case II, is discussed here for two reasons: it illustrates 

the effect of structural elasticity on wheel shimmy, and it is better 
suited than case I for evaluating approximations Dl, D2, and D3 as applied 
to landing-gear problems involving structural elasticity, (it may be 
recalled that these three approximations were of little value in dealing 
with a rigid landing-gear strut (case I) ; however, for a flexible strut 
these approximations may sometimes be of value.) In discussing case II, 
and also later cases, no further mention will be made of the previously 
published theories or of the question of agreement between theory and 
experiment; all discussion will be restricted to the summary theory and 
its systematic approximations. 

The discussion of case II proceeds as follows: The equations of 

motion for case II are derived according to the summary theory And its 
systematic approximations. As for case I, it is more convenient to rede- 
rive these equations of motion in a slightly different manner rather than 
to apply the equations derived earlier for the completely general case. 

After these derivations are made the equations for the stability bound- 
aries are established. Finally, the damping required to prevent shimmy 
is presented in curves as a function of strut stiffness and rolling velocity 
for a specific sample landing-gear configuration according to the predic- 
tions of approximations C, Dl, D2, and D3. (For the present case approxi- 
mations Cl and C2 are identical and are henceforth referred to collectively 
as approximation C.) These curves are utilized to obtain some insight into 
the accuracies of the predictions of approximations Dl, D2, and D3 with 
respect to the more, advanced approximation C. 


General Derivation 

The derivation of the equation of motion for the summary theory pro- 
ceeds as follows. The details of the landing gear considered are illus- 
trated in figure 6. This gear has a rigid symmetrical swiveling part 
having a mass m and a moment of inertia about its center of gravity I Q . 
The nonswiveling part of the landing gear consists of a spring of stiff- 
ness k with an attached mass m^. The lateral displacement of the 
swivel axis is designated as n_ . 

'cl 
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Setting the sum of the lateral spring and inertia forces acting on 
the swiveling part equal to the inertia reaction of its center of gravity 

mDt 2 ( T )a “ c 2® ) yields the relation 

Vo “ Hi “ m l D t\ = m D t\ " mc 2 D t 20 Vo) 

Substitution for Aq from the relation 

^0 = yo - = y 0 - ^a + ae 

(see fig. 6) yields, after rearrangement, 

Vo " ( m l D t 2 + m D t 2 + \ + k ) T1 a + ( mc 2 D t 2 + aK A) 0 = 

Setting the sum of the moments about the center of gravity 
swiveling part equal to the inertia reaction yields the result 

K a a " ¥o C l “ kT k c 2 “ m l D t 2l k c 2 - g D t 0 - P0 - tv D t A Q = I Q D t 2 9 (153) 

(see fig. 6) where I 0 represents the moment of inertia of the swiveling 
structure at its center of gravity (l Q = 1^ - mc^ 2 ). Substituting for 

a and Aq according to equations (23) and (151) then yields, after 
rearrangement , 

(tv D t - K a v“\ + K ACl )y 0 + (m^ D t 2 - rv D t + kc 2 - K^cJt^ + 

^I Q D^ 2 + g D t + Tav D t + P + K a + ac-jK^Q = 0 ( 15^ ) 

The third equation for this system for the general case is given by 

the kinematic relation of equation ( 19a) or (19b). When 7 is omitted, 

the space derivatives are replaced by time derivatives, and tjq is set 

equal to q a - a0, this relation becomes 

-^1 + z-jv -1 D. t + z 2 v ^t 2 + * • -)y 0 + (^1 - a )s + Tl a = 0 (155a) 

from equation (l9a) or 

-(l + Lv _1 D t )e hV Dt y Q + (l 1 - a)0 + T] a - 0 (155b) 

from e quat ion ( 19b ) . 


(151) 


0 (152) 

of the 
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The three equations (152), (154), and (155&) or ( 155b ) completely 
describe the motion of the landing gear according to the summary theory 
in terms of the three variables Yq, T) a , and 0. The corresponding equa- 
tions for the systematic approximations are obtained as follows. 


Systematic Approximations 

Approximation A .- For case II, the three governing equations of 
motion for approximation A are the force and moment equations (152) 
and (154) and the kinematic equation: 

-(1 + IqV - ^ + l 2 v' 2 D t 2 + l 5 v" 5 D t 5 )y 0 + - a )e + Tl a = 0 (156) 

Equation (156) is obtained by omitting all l n 's T° r n> 3 in the gen- 
eral kinematic equation (l55a). 

Approximation B .- The three governing equations of motion for approxi- 
mation B are the force and moment equations (152) and (154) and the kine- 
matic relation 


-(l + Zqv" 1 !^ + Z 2 v _2 D t 2 )y 0 + [l ± - a)e + q a = 0 (157) 

which is obtained by omitting the l ^ term in equation (156) . 


Approximation C .- The kinematic equation for approximation C is 
obtained by omitting the lg term in equation (157) • The resulting 

relation is 


-(l + V _1 D t )y 0 + (z-L - a) 


6 + q = 0 


(158a) 


or 


^a 





(158b) 


The force and moment equations (152) and (154) also apply for this approxi- 
mation. However, a slightly simpler form of these equations can be obtained 
by substituting for q a , according to equation (158b), in the terms con- 
taining K A T] a and ttj^ in these two equations. With this substitution, 
the force equation becomes 
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-Nv^^D^yQ - + m D^ 2 + + ^mc 2 D^. 2 + N^9 = 0 

and (using the relation = eN) the moment equation becomes 


(159) 


N-TDt 

K A 


- Nv - ^(c-j_ + e^D.j- 




+ ( mi c 2 D t 2 + kc 2 )q a + 


p NtvD + / \ 

I o D t + S D t + — + p + (cq + e ) N 


6=0 


(160) 


Equations (158a ) } (159) ( or (152) )j and (160) (or (154)) describe the 
motion for approximation C. 

Approximation Dl .- The equations of motion for approximation D1 are 
obtained by setting ~Kq equal to 0, or 


and 


yo = % " a0 


K A »« 


(161) 


in the force and moment equations (159) and (l60) . Thus 

-^m-^ D^ 2 + m D^ 2 + Nv *4)^ + kjq a + ^mc^ D^ 2 + aWv _1 D t + N^0 = 0 (l 62 ) 


m l C 2 D t 2 " Nv_1 ( c i + e ) D t + kC £ 


aNv” 1 (c n + 

nD, + p + (c., + 

An 

\ 1 

/ t K V 1 

J 


K + 


0 = 0 


X o D t + 8 D t + 


(lb3) 


Equations (162) and (163) completely describe the motion for approxima- 
tion Dl. 

Approximation D2 . - The equations for approximation D2 are obtained 
as follows. In the force equation (152) set a equal to 0, or 


0 = v -J "D, y 


t J 0 


(164) 


which gives the relation 
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^mc 2 v -1 D t 5 + aK^v“ 1 D t + K^y Q - D t 2 + m D t 2 .+ + k^T} a = 0 (165) 

For the moment equation ( 153 ), set a equal to -Kq (eq. (85)) 

according to approximation C, and apply equations ( 151 ) and (l 64 ). The 
result is 


I 0 v + gv" 1 D t 2 + aT D^ 2 + tv + pv -1 !^ + a^c-^ + e^I^v -1 ^ + 


( C 1 + 6 K 


y o + 


m l C 2 D t “ TV 


D t - ( C 1 + e K + c 2 k ] T l a = 0 


(166) 


Equations (165) and (l66) are the basic equations for approximation D 2 . 


Approximation D 3 «- The basic equation for approximation D 3 is obtained 
by first solving the equations (162) and (163) of approximation D 1 simul- 
taneously to eliminate either rj a or 0 and then letting N approach 00 

so that all terms not multiplied by N vanish. The resulting equation, 
after dividing out the factor N and using the relation c q + c 2 = a > 

can be expressed in terms of 0 as 


^ I Q + m^a(a + e) + mc^(c 1 + e)J D_^ + - |m^(a + e) + 

m(cp + e)Jv + g >D^ 2 + jp + ak(a + e)J + vk(a + e)je = 0 ( 167 ) 


Stability Boundaries 

The stability boundaries for case II are obtained in the same manner 
as those for case I. For the summary theory they are obtained as follows. 

Purely oscillatory boundaries .- The equations for the purely oscil- 
latory motion boundaries are obtained by substituting into the differen- 
tial equations the expressions 


0=0 


max 


ivt 


"a ‘ V 


i(vt+cn) ivt. 

- L/ = T) e (cos pj. + 1 sia °i) 
jiiax 


y 0 = y (W 


i(vt+a2) ivt 


x 


= yg ^e (cos + i sin 02) (l68) 
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Substitution of these relations into equations (-152), (15 ^) , and ( 155b ) , 

differentiation and cancellation of e : '" yt and separation of real and 
imaginary parts into separate equations yields the expressions: 


From equation (.152) , 


K A (y 0 cos 
A \ u max 


a 2 ) + (m-jV 
( aK A - me 2 v 2 )e max « 0 


2 + mv 2 - K-v 


-)( 


Snax ° 0S °1 + 


(169) 


and 


M yo max Sl " ° 2 ) * (“l'' 2 + mv2 ' ^ " k X%K Sin °l) = 0 (lT0) 


From equation (15 1 +)^ 
Ci K- 


“ s ° 2 ) - ( tvv - sin ° 2 ) 


(- 


m i G 2 v ^ + c 2^ 


- °A)(V COG °l) + TVV ( T Sax Sln °l) + 


(-I 0 v 2 + P + K a + aCjKjj 


e =o 

max 


(171) 


and 


c A(»V Sln ° 2 ) * ‘ V'MKbx °° S S s) + ( 


-mc2V + 


k0 2 ' Vl)(v Sln °l) ‘ TVV (%ax °° B °l) + ( gv + Tavv ) 0 Jnax = 0 


(172) 


From e quat ion ( 155b ) , 


-P 


lco 


Kax C0S “ 2 ) + P2 “ ( y Vx Sln ff a) + (h - a ) 0 max + (%« =° 5 «i) - 


0 


( 173 ) 
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and 


- p 2»Kax C ° S °2) - p l»( y 0, 


max 


sin 0 \ + T) a sin 0=0 (17I4-) 

2/ Tnax 1 


Equations (169) to (17*0 can "be considered as six linear sim ult aneous 
algebraic equations with no constant terms in the five variables 


y 0 cos y 0 

u max 2-' o, 


sin 0 


max 


2 J T * a max 


COS 0 


1 > 


sin 0]_, and 8 


max 


Then for this system of equations to have solutions other than zero it is 
necessary that the determinant of the coefficients of any group of five 
of these six equations should equal zero. The determinant for equa- 
tions (169), (170), (171), (173), and (174) is 


'A 

k a 

0 

-Pi 

-p 2 


-1 2 

-tvv + K a v v ~ m i c 2 v + c 2 k - 


h. 

p 2 

-Pi 




v - K-x - k 


(“i + m ) 


0 

1 


-V + P + ^ + ac l K A 


" 2^ 
0 

2-i “ & 


. 0 (175) 


and the determinant for equations (169), (170), (172), (173), and (174) 
is 


■vrv - y -1 v c 1 K x 

k a 0 

0 k a 

-pi 52 

-p 2 -Pi 


-rvv 

+ m)v 2 - - k 

0 
1 
0 


2 

-hk^v + keg - ^\ c i 
0 

(m^ + m)v 2 - - k 

0 
1 


gv + raw 
aK^ - mcgV 2 
0 



0 


= 0 (176) 


where = p loo and p^ = p^ for the summary theory. The corresponding 

equations for approximation A are obtained by setting p equal to 1 - l 

3 1 

and p^ equal to ~ > for a PP roxima tion B by setting p^ 

2 

equal to 1 - Vi and Pg equal to Vi< and for approximation C by 
setting p-^ equal to 1 and pg equal to. 

The two equations (175) and (176) completely describe the conditions 
governing purely oscillatory modes of motion according to the summary 
theory and approximations A to C (other groupings by fives of the equa- 
tions merely lead to repetition of these two relations). The procedure 
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for obtaining the stability boundaries for each of the other systematic 
approximations (Dl, D2, and D3) is similar to that just outlined. The 
resulting stability-boundary equations for the other approximations are 
listed in appendix E. 

Purely uniform motion .- For purely uniform motion all variables will 
have constant values which may be represented as 


T * a max 

y o = y o mx 

Substitution of these relations into equations (152), (15 ^) , and (155) 
yields the results 


aK *Vx ' ( K X + k )%ax + VCW = 0 
(p + + “l K A) f; ,r,ax + ( c 2 k - c A)Vx + 'A’V ' ° 

(h - a)0 max + V ' J Vx = ° 

For nonzero solutions of these three equations the determinant of the 

coefficients of 0 , tu , and y n must be zero. Evaluation of 

max- 7 Ia max u max 

this determinant gives simply 

a + e + | = 0 (177) 


Character of the Motion Between Stability Boundaries 

In order to determine the character of the landing-gear motion (sta- 
ble or unstable) between stability boundaries it is first convenient to 
solve the equations of motion for each approximation simultaneously to 
obtain a single linear differential equation in one variable for each 
approximation. From these differential equations the stability of the 
motion may be determined by examining the corresponding characteristic 
equations by any of the methods discussed in appendix C. These charac- 
teristic equations for case II, according to the various systematic 
approximations, are listed in appendix F. 
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Evaluation of Approximations Dl, D2, and D3 

In the earlier discussion of case I it was not possible to present 
a fair relative evaluation of the three parallel approximate theories Dl, 
D2, and D3 since for case I none of these theories provides any realistic 
information. However, for case II such a comparison can he made between 
the predictions of these three approximations and the more accurate 
approximation C, and a specific example will be discussed here for a sam- 
ple landing-gear configuration having the relative dimensions and proper- 
ties L = 0.8r, h = a = 0.5r, c p = c 2 = 0.25r, 6 = 0.3r, m-j_ = 0.35m, 

I Q = mr^, and t = p =0. The actual calculated behavior of this landing 
gear in terms of damping required for stability as a function of rolling 
velocity according to approximation C is shown in figure 19 for four 
values of the ratio of strut stiffness to tire stiffness k/K^. It is 

seen from this figure that as the stiffness of the strut is decreased 
from infinity the damping requirement is increased. Also, for large 
values of strut stiffness the region of maximum damping required lies at 
low speeds, whereas for small values of strut stiffness it lies at higher 
speeds . 

The theoretic predictions of the three theories Dl, D2, and D3 
for this sample landing gear are compared in figure 20 with the corre- 
sponding predictions of the more accurate approximation C (from fig. 19) 
for three values of strut stiffness, k = 0.2K^, 1.0K , and 5*0K^. It 

is seen that for each strut stiffness approximations D2 and D3 provide a 
considerable overestimate of the damping required for stability. On the 
other hand, approximation Dl gives results in good agreement with those 
of approximation C for the ratios k =0.2 and 1.0, but this approxi- 
mation greatly underestimates the damping for the large value of strut 
stiffness k/K^ = 5 •Cl- 
in view of the comparisons of figure 20, it appears that approxima- 
tions D2 and D3 will not, in general, give reliable quantitative estimates 
of the damping required for stability. It appears that approximation Dl 
may give reasonable results for some cases in which the lateral stiffness 
of the strut does not greatly exceed the lateral stiffness of the tire. 
Since this latter conclusion is based on only one set of landing-gear 
parameters, the degree to which it is valid in general would require a 
more extensive investigation for a range of landing-gear properties. 


Practical Application 

One limitation on the practical application of the preceding equa- 
tions for case II lies in the assumption that the damping is simply pro- 
portional to the angular swiveling velocity D^.0. As was previously 

mentioned, Moreland has demonstrated that this assumption is sometimes 
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unreliable since it implies a neglect of the torsional flexibility of the 
landing-gear strut, which in turn can sometimes lead to a false prediction 
of stability for heavily damped systems (see ref. 11 or 12). Thus, for 
systems in which torsional flexibility of the strut is important, it will 
be necessary to replace the damper unit of case II by a damper and spring 
in series, as has been done by Moreland (refs. 11 and 12), where the 
spring represents the strut torsional stiffness. This particular case 
of a series damper-spring unit applied to the landing gear of case II is 
not considered separately in this paper; it is, however, included in the 
more general case III to be considered next. 


CASE III 


The next type of landing-gear construction to be considered is chosen 
largely to illustrate the application of the summary and approximate 
theories to more complex problems than have previously been considered 
by now making use of transfer-function concepts. This landing gear is 
assumed to be of the same general type as that of case II except for the 
following generalizations. In case II the lateral deflection character- 
istics of the landing-gear strut were represented by a single spring and 
mass combination, or, more precisely, the force exerted on the swiveling 
part of the landing gear by the strut was set equal to 

F n = -(k + mi D t 2 )T| a (178) 

For case III it is assumed that the strut (or, more generally, the sup- 
porting structure) is a more complex linear system than is a spring-mass 
combination, so that the strut force-deflection relation of equation (178) 
can be generalized to the new form 

= -hoKK hra) 

where T-^q is a function of the differential operator which repre- 

sents the transfer function correlating F^ and T] a *^ 

Similarly, the moment on the swiveling part due to the strut, which was 
previously set equal to 


Mg = -( P + g Dfc)0 


( 180 ) 


^Moreland, in reference 11, has advanced a similar generalization of 
the strut lateral-deflection characteristics by means of a concept of 
"virtual elasticity." However, Moreland’s generalization is less general 
in that it does not provide for the existence of strut structural-damping 
forces . 
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for case II, is now generalized to the form 

% = - T ii( D t) e ( l81 ) 

The differential equations for case III are easily obtained from the 
corresponding differential equations for case II by replacing k + m^ D t 2 
by T 10 an d p + g by T - 1 n . For example, for approximation D1 equa- 
tions (162) and (163). are replaced by the new relations 


m D^ 2 + Nv""^D-j- + T-j_q 


( D t)] ^a + (’ 


me g D t 2 + aNv~ 1 D t + N]0 = 0 


je = 0 (182) 


+ e) Dt + c 2 T 10 (D t ) q a + I Q D t 2 + aNv" 1 ^ + e)^ + 


( C 1 + e ) 


N + T 


11 


( D t) 


0=0 


(183) 


which lead to a characteristic equation of an order depending on the order 
of the T's. If the T's are analytically defined functions, the calcu- 
lation of the stability of the motion and the stability boundaries pro- 
ceeds exactly as for case II. However, if the T's are not analytically 
defined functions (for example, if they are determined by experimental 
tests) a slightly different procedure of the following type is needed. 


In order to determine partial information about T-^q, the swiveling 
part of the landing gear can be removed and the remaining strut can be 

ivt 


subjected to a periodical lateral force 


•^rj _ ^ti e 
1 max 


(either by cal- 


culation or by actual vibration tests) . The resulting lateral-deflection 
response of the structure will have a certain amplitude and phase shift 
which are given by the relation 


•n = — F 

T 10 (iv) 1 


(184) 


which is obtained by substitution of the sinusoidal variation for F^ 
into equation (179) • The function 1/T^g(iv) is a complex function of 

the circular frequency v, the absolute value of which represents the 
amplitude response and the argument of which represents the phase shift j 
it is generally called the frequency-response function of the system. 

Similarly, a f re quency-re sponse function is defined for the response 
of the landing-gear strut to torsional moment oscillations by the relation 
(see eq. (l8l)) 



98 


MCA TN 3632 


Q = - m % (185) 

Tll(xv) 

With the aid of the experimental or calculated functions T-^q(Iv) 
and Tn(iv) , the stability boundaries for any of the theories may be 
obtained by the usual procedure of substituting expressions of the form 

e iy k into the corresponding differential equation together with the T's. 
For example, for approximation D3 the basic differential equation is 


I o + mc l( c l + e )] D t 3 + 

v(a + s)T 10 (iv)| 0=0 


m( c i + e)vjD t 2 + l^ppCiv) + a(a + e)T 1Q ( iv) + 

( 186 ) 


(obtained by converting eq. (167) to apply to case III), and the stability 
boundaries for purely oscillatory motion (obtained by setting 9 equal 

to e ivt in equation (l86) and separating real and imaginary parts) are 
given by the simultaneous equations 


-m^c 1 + e)vv 2 - vl T-y^iv) + a(a + e)T 10 (iv)J + v(a + e)R T 10 (iv)J = 0 

(187) 


-I Q v - me 


v 


1 (c 1 + + vR jl 11 (i v ) + a ( a + e ) T io( iv |] + 

(a + e)l|T 10 (iv)J = 0 


(188) 


where R and I represent the real and imaginary parts of the bracketed 
functions. Analogous equations are obtained in a similar manner for the 
other approximations. 


In regard to the question as to whether any particular motion between 
stability boundaries is stable or not, case III may present a more diffi- 
cult problem if the forms of the T-functions are not known in terms of 
ratios of polynomials, that is, if only the frequency-response variations 
are known. In this event, for example, the usual form of the Routh-Hurwitz 
stability criteria (which is applicable to polynomial forms only) cannot 
be applied and criteria of the Nyquist type must be used. A brief dis- 
cussion of these criteria is contained in appendix C. 

The procedure for applying the summary and systematic approximation 
theories to cases of arbitrary complexity is essentially the same as the 
procedure discussed above for case III, the only important difference 
being that for the general case the equations of motion (eqs. (l6) or 
(19), (63), (6k ) , (65), (67), (69), and (71)) are more numerous and more 
complicated. No new concepts need to be discussed. 
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CONCLUDING REMARKS 


This paper has presented a correlation and evaluation of the pre- 
viously published theories of linearized tire motion and wheel shimmy 
and has demonstrated that the major merits of all of these theories are 
contained in a summary theory which represents a minor modification of 
the basic theory of Von Schlippe and Dietrich. In cases where there are 
strong differences between the existing theories and the summary theory, 
the previously published theories have, in the main, been demonstrated 
to possess certain deficiencies except for Moreland* s advanced theory, 
for which no adequate evaluation was possible. 

A series of systematic approximations to the summary theory has been 
developed herein for the treatment of problems too simple to merit the 
use of the complete summary theory. These systematic approximations have 
been shown to resemble closely the previously published theories except 
that in some details they avoid some of the limitations encountered in 
these theories. 

Comparison of the existing experimental data with the predictions 
of the summary theory and its systematic approximations has indicated a 
fairly good degree of correlation between the higher approximations and 
the existing experimental data for the cases investigated. However, 
since the agreement is far from perfect in some respects and since most 
of the limited amount of existing experimental data was obtained with 
small models there still remains the question as to whether the theory 
is safely applicable to full-scale conditions. In particular, the impor- 
tance of hysteresis damping remains undetermined. 

In regard to the determination- of the various tire constants required 
for theoretical shimmy calculations, it is noted that the existing perti- 
nent experimental data, mostly contained in references 21, 2 h, and 32 
to 37, are extremely limited and apply mostly to small, obsolete, or 
foreign tires. Furthermore, although various attempts have been made 
to correlate and to reconcile theoretically the experimental data (e.g., 
ref, 2 or 38 ), there still apparently does not exist any fully reliable 
theoretical means for predicting all of the needed elastic characteristics 
of tires. In view of these considerations, a need exists for additional 
experimental data on modern tire characteristics and also for a more ade- 
quate evaluation of the existing data to determine whether these data can 
be applied by scale laws to predict the characteristics of any tire with 
tolerable accuracy. 

In regard to the adequacy of a linearized theory of tire motion, it 
can be stated only that there is as yet no strong indication that a non- 
linear theory is required for prediction of the stability boundaries. 

If, however, a knowledge of the large-angle (nonlinear) behavior is 
required, a theoretical system for dealing with this problem could be 



100 


NACA TN 3632 


developed on the basis of assumptions of the type advanced by Kelley, 
Rotta, and Temple (see refs. 18, 19, and 21 (p. 36 ), respectively). 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., January 13, 1956. 
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APPENDIX A 


CALCULATION OF EQUIVALENT VISCOUS DAMPERS 


GENERAL CASE 


In the derivations and equations in this paper, only linear damping 
ter ms were introduced so that the resulting equations would remain lin- 
earized. However, the damping moments caused hy friction,, hysteresis 
losses, and the ordinary shimmy dampers are nonlinear and therefore it 
is necessary to replace these nonlinear moments with equivalent linear 
viscous moments. The equivalent viscous moment is usually determined 
by ass umi ng that linear and nonlinear damping moments are equivalent if 
they dissipate the same amount of energy during each cycle of shimmy 
oscillation. 

For a linear damper of moment g D-^ijr the energy dissipated per cycle 
of sinusoidal oscillation is (for \|r = T| f mav sin vt) 

r> 2rt 

E = J g D t t di|r 


and, since D, i|r = t | f v cos vt and di|r = l ay v cos vt dt, 

' t max max 

>2n 


0 

E = gv\lr / cos^vt d(vt) 

B Y max J 0 


which gives 

E = ^S^max 2 


Therefore the linear damping constant g is related to the energy loss 
per cycle by the relation 


g = 




max 


(Al) 


By using this relation an effective value of linear damping constant can 
be calculated for any nonlinear damper if the energy dissipation per cycle 
is known. 
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VELOCITY-DEPEKDENT DAMPING 


By using equation (Al) Rotta (ref. 2) has calculated the effective 
damping constant for damping moments of the type 


= l M o + 


n l\ D t t 

t Y l /IeT 


(A2) 


The first term represents friction damping and the second represents 
fluid damping. The exponent n^ will probably always be between 1 and 2. 

Rotta' s calculations proceed essentially as follows: First, the 

total energy dissipated is calculated from the relation 


+1 


E 


max 


0 




(A3) 


where 


t = W sin vt 


(A4) 


After combining equations (Al) to (A4) and integrating, the following 
equation for g is obtained: 


g 


n-, -1 n-i-1 k P rt / 2 n-,+1 . . 

— 2 — + Mmax 1 v 1 ^ / cos 1 vt d(vt) (A5) 


jtvt-. 


max 




0 


It is seen that the damping constant depends on the amplitude 

for n^ > 1 it is large at small angles due to the first term and is 

large at large angles due to the second term. The minimum value of the 


/ dg 

equivalent damping constant (obtained by setting 

\ d ^max 

occurs at the angle 

M, 

%. = 


equal to 0 


n-i / \ o ^ /2 


n-,+1 

cos vt d(vt) 


l/ni 


(A6) 


and the corresponding minimum damping constant becomes 
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where 


g 


n -^-1 



l/ n l 


k n = 


4n n 


ni-1 

rt ( n l " x ) ni 


,:It / 2 nq+1 . 

cos vt d(vt) 


0 


l/n-j 


(A7) 


(A8) 


For the special case of velocity- squared damping (n = 2) equation (A7) 
reduces to the relation 



(A9) 


FRICTION DAMPING 


The equivalent damping constant for friction damping (constant 
moment M 0 ), obtained by setting q n equal to 0 in equation (A5 ) , is 


g = 


4M„ 


itvt 


max 


(A10) 
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APPENDIX B 

DIFFERENTIAL EQUATIONS FOR CASE I WITH INCLUSION 
OF HYSTERESIS EFFECTS 


If the differential equations of approximations B and C2 for case I 
(see eqs. (122) and (124 ) , respectively) are modified to take into account 
the hysteresis force and moment expressions derived in this paper, the 
following differential equations are obtained: 


For approximation B, 

E o d1 % + E i d3 yo + E 2 d2 Yo + e 3 Dy o + E 4 y 0 = 0 


where 


(Bla) 


E 0 = V 2i 2 

E-j_ = I 1 | f v 2 Z 1 + (a-rv 2 cos k + gv)lg + (a 2 K^T^ + cos^ZgV 

Eg = Iv 2 + (arv 2 cos k + gv)i- L + (a 2 K^ + cos 2 k -k p + P K )^ 2 + 

(a 2 NT-> + a'K T_ cos k) v 
\ A a a ) 


E^ = a 2 N + a'K a cos k + pZ^ + P^Z-^ + gv + 0 Z-^tv 2 cos 2 k + 


(craNT-,. cos k + K T cos 2 k) v 
\ A a a / 

Ek = craN cos k + K cos 2 k + p + p„. + u„ 

CL ft lx 


> (Bib) 


For approximation C2, 
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where 

E o = V 2l i 

E 1 = I f v2 + ( aTv2cOS K + gv)li + (a 2 KT-^ + ai ^(, T a cos K^V 
E 2 = a 2 N + aK^ cos k + pZ^ + P K ^]_ + gv + Z^tv^cos^k + 
(aNT^ cos k + I^T a cos^j v 
Ej = aN cos k + cos^k + p + p K + 


!> (B2b) 




io6 


NACA TN 3632 


APPENDIX C 


STABILITY CRITERIA 


In this appendix a brief review is presented of most of the existing 
methods for examining the stability of motion of systems whose motions are 
governed by linear differential equations of the type 

f(Dt)9 = 0 (Cl) 

The solution of this type of equation consists of terms of the form 

8 = e pt (C2) 

whence 

D t n 0 = p n e pt = p n 0 (C3) 

Substitution of equation (C3) into equation (Cl) yields the algebraic 
equation 

f(p) = 0 (c4) 

for the p's. Equation (C4) is called the characteristic equation of 
the differential equation (Cl) . 

The type of motion for the linear system is determined by the char- 
acter of the complex roots of the characteristic equation. Most impor- 
tant, the motion is entirely stable if and only if the characteristic 
equation possesses no roots having positive real parts. Several proce- 
dures are available for determining whether a particular characteristic 
equation has such roots with positive real parts. 

One procedure which is useful in cases where the characteristic 
equation (C4) can be written in the polynomial form 

a Q p n + a-jp 21 " 1 + . . . + a n = 0 (C5) 

is the well-known Routh-Hurwitz criterion which makes use either of the 
Routh test functions (ref. 26) or of the equivalent Hurwitz determi na nts. 
In Hurwitz' form the requirement for stability (or no roots with positive 
real parts) is that a Q and all of the n determinants 

D 1 = a ! 
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. 


a, 


2u-l 


a 0 a 2 a 4 ’ ' ' a 2u-2 


D u - 


0 3 * 2 ^ • • * ^2u.°“3 


® a 


• « • 3 -, 


U 


for u = 1, 2, . . . n must be greater than zero, or, in the alternate 
form of Cremer (ref. 39) y all of the a n 's and either the even or -the 
odd Hurwitz determinants must be positive. 


This criterion is particularly suited for examining the stability 
of linear systems with polynomial-type differential equations of low 
order. However, for high-order polynomial -type differential equations 
this criterion may not be the easiest to use and for nonpolynomial-type 
equations the criterion is not directly applicable. For such cases use 
can be made of the graphical-type criterion originated by Nyqu'ist 
(refs. 28 to 31)- Some discussion of criteria of this type is contained 
in most books dealing with servomechanisms or feedback amplifiers (for 
example, refs. 29 and 31)- These references provide the necessary theo- 
retical information for applying these criteria and the theory will not 
be repeated herej however, it may be useful to set down here, together 
with an example, one mechanical procedure for applying this criterion 
and a few pertinent comments. 

Consider a differential equation with the characteristic equation (b 4) 
for a case where the function f(p) cannot necessarily be easily expressed 
in a simple polynomial form which can be handled by the usual Routh-Hurwitz 
criterion. (This may be the case, for example, where part of the func- 
tion f(p) is evaluated from experimental frequency-response data.) 

The function f(p) is assumed to be a single-valued function of p 
which is real when p is real. It is also assumed that the function f(p) 
has no poles in the region of the complex p-plane where the real part of 
p is greater than zero. When the equations of motion are set up in the 
manner followed in this paper, the condition of no poles in this region 
is usually satisfied for actual landing gears since this condition implies 
only that the landing-gear strut, as represented by equations or experi- 
mental curves, possesses some damping or is at least not inherently an 
unstable structure. For example, for the equations of case III, only the 
poles of the functions T 1Q (p) and T-q(p) could lead to such poles. 

However, if, for example, T 1Q (p) had such a pole, equation (l6l) would 
indicate the possibility of a steady or divergent oscillating force F Tj 
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corresponding to this pole, even if the lover. end of the landing-gear 
strut vere held fixed (q a = 0)j this obviously cannot occur in actuality. 

In order to decide whether the motion of the system described by 
equation (C4) is stable, the following procedure may be followed: 

(1) Determine the variation of f(p) in equation (C4) for the case 
of pure sinusoidal oscillations, that is, for p = iv in the ra ng e 

0 < v < 

(2) Plot the real part of f(iv) against the imaginary part for 
the complete range 0 < v < <». This will give a curve such as is illus- 
trated in figure 21 for a sample case. As v varies from 0 to °° 
this curve will move about the origin through a net angle of Jrt radians, 
this angle being considered positive in a counterclockwise sense, (j = 1 
for the case illustrated in fig. 21.) 

(3) Determine the asymptotic behavior of the characteristic func- 
tion f(p) or f(iv) for p °° or v— at this limit the function 

will behave as f(p) oc pJ or f(iv) « vJ, whence j can be determined. 

(4) Under the preceding restrictions f(p) being a single-valued 
function of p, real when p is real, and having no poles in the half- 
plane R(p) > 0) , it can be shown that the motion corresponding to the 
differential equation is stable if and only if 2J = j. (in the sample 
case of fig. 21, where J = 1 and j = 2, the motion is therefore 
stable . ) 

In conclusion it might be noted that, although either the preceding 
Nyquist type stability criterion or the Routh-Hurvitz criterion can 
usually be applied to most of the approximate equations discussed in 
this paper, they cannot be directly ■ applied to some equations of trans- 
cendental form such as those of the summary theory, since such equations 
may correspond to infinite-order linear differential equations (for 
example, see eq. (120) ) ; hence, an infinite number of Hurwitz determinants 
would have to be evaluated or the Nyquist type plot would circle the 
origin an infinite number of times. 
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APPENDIX D 

STABILITY BOUNDARIES FOR CASE I 

The following equations describe the conditions under which purely 
oscillatory motion is possible for case I for the summary theory and the 
systematic approximations. 

For the summary theory and approximations A to. C2, 



where, for the summary theory, 


Pl = p loo = cos Vqh - Lvj_ sin Vqh 
Pg = P 2 ro = sin Vqh + Lv-q cos Vqh 

for approximation A, 

i 7 2 

P 1 “ 1 " Z 2 V 1 

p 2 = l l V l " Z 3 V 1 5 

for approximation B, 

7 7 2 

Pq — 1 “ ^2 V 1 


p 2 ■ Z l v l 
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and for approximations Cl and C2, 


P 


1 


= 1 


= Vi 


For approximations D1 and D3 purely oscillatory motion does not 
exist. 


For approximation D2, 



a 2 K^ + aeK^ cos k + p + p K + tv^cos 2 k 



(D3) 


g = v 


1^ (aR^ cos k + eK^ cos 2 k + c^F z sin k cos k) 
a^K^ + aeK^ cos k + p + p K + tv 2 cos^k 


ar cos k 


m 


The stability boundaries for -uniform motion are obtained by setting 
the coefficient of the yQ terms in the various differential equations 

equal to zero. For example, for the summary theory and approximations A 
and B the equation 

O 

oaW cos k + cos k + p + p K + u K = 0 


describes this stability boundary. 
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APPENDIX E 

STABILITY BOUNDARIES FOR CARE II 


The equations governing purely oscillatory modes of oscillation for 
case II are as follows for approximation Dl, D2, and D3: 


For approximation Dl, 


V s = v 2 !) 2 


Jr o + mc^c-^ + c) + m-^afa + € ) j a 2 - jp + ak(a + ejj 


[(mjL + m)v 2 - kj j"- |l 0 ( ra Q +■ m) + mm 1 c 2 2 Jv l4 ‘ + jl Q k + mkc 2 2 a m-jN{a + e) + + e) + + m)jv 2 - kjp + N(a + e)] 


(El) 


g = - Nv 


-1 + + m^(a + e )J - k( a + e)^ Jj^m^a + mc p) v 2 + N - akj 

j^m^ + m^v^-k + N^v" 


-2 2 
r c -y c — 


(E2) 


For approximation D2, 


jl Q (m 1 + m) + mm-jCg 2 ' V* - |l 0 (K^ + k) + p(mi + m) + mc 2 2 k + mjK-^afa + e) + mK ?v c 1 ^c 1 + e)Jv 2 + |p ( k^ + k) + kK^a(a + e)J 


TV 2 (m^ + m) - Tk 


(EJa) 


-V -i 


( Tam, + Tc-,m 


1. 'T 


) - jaTk + m-jJK^ ( a + e ) + mK^ (c ^ + e)Jv 2 + kK^(a+ e) 1 


For approximation D3, 


v 2 = 


v 2 r(rrii + m ) v2 " ( k a +k )J 


p + ak(a + e) 


g = v 


I Q + m^a(a + e) + mc q( c i + 
k(a + e)v ^ - m^c^ + e) - m^(a + e) 


(E3b) 

(E4a) 

(E4b) 
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APPENDIX F 

CHARACTERISTIC EQUATIONS FOR CASE II 

The characteristic equations for case II are as follows for the 
summary and systematic approximation theories: 

For the summary theory and approximations A to C2, 


k A 


-( m l + “)p 2 - K A - k 


mc 2 p 


+ a-K-v 


( tv - K jx v” 1 )p + CjK^ m^gp 2 - -rvp + ke 2 - K^c 1 IqP 2 + (g + Tav)p + p + + ac-jK^ 


= 0 


-f(p) 


1 


l l 


a 


where, for the general theory, 

f(p) = (l + lp 

for approximation A, 

f(p) = 1 + i 1 v -1 p + i 2 v _2 p 2 + ZjV 5 p3 

for approximation B, 

f(p) = 1 + + ^ 2 V ” 2 P 2 

and for approximation Cl or C2, 


f(p) = 1 + i 1 v“' 1 p 
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For approximation D1 the characteristic equation is 

^(mp + m) + xnm 1 c 2 2 Jp it ’ + jg(mp + m) + Nv -1 jl 0 + m-ja(a + e) + 

mc l( c l + e )J > + [(lo + mc 2 2 )k + p(m 1 + m) + mN^cp + e) + 

m-jll(h + e) + gNv"-^Jp 2 + <fja(a + e ) Nv_1 + sj k + pNv” 1 
|p + (a + e)Njk = 0 


For approximation D2, 


|l 0 (mp + m) + + j^rm^a + me-^ + g^m- jL + m)v ^ 


P 4 + 


jl 0 ^K^ + k^v -1 + p(mp + m)v -1 + rv^mp + m^ + kmc^v -1 + m^a + e)K^v -1 + 
mc 1 (c + e)K^v”'*" p + jjak + g(K^ + k^v -1 + m^ + e)K^ + m-^a + e)K-^Jp 2 + 

jjrvk + p(% + k)v _1 + akK^(a + e)v~^jp + kK^(a + e) = 0 


For approximation D3 see equation (167). 
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TABLE I 


PRIMARY ASSUMPTIONS FOR THE VARIOUS THEORIES OF TIRE MOTION 1 


Theory 

N 

K X 

Ka 

e 

*1 

l z 

Z 3 

M n > 3) 

Remarks 

Summary theory 

F 

F 

F 

F 

F 

F 

F 

F 


Approximation A 

F 

F 

F 

F 

F 

F 

F 

0 


Approximation B 

F 

F 

F 

F 

F 

F 

0 

0 


Approximations Cl and C2 

F 

F 

F 

F 

F 

0 

0 

0 

- 

Approximation D1 

F 

00 

F 

F 

0 

0 

0 

0 

Assumes laterally rigid tire 

Approximation D2 

00 

F 

00 

F 

00 

0 

0 

0 

Assumes torsionally rigid tire 

Approximation D3 

00 

00 

00 

F 

F 

0 

0 

0 

Assumes laterally and torsionally rigid tire 

Von Schlippe-Dietrich and Rotta 

F 

F 

F 

F 

F 

F 

F 

F 


Bourcier de Carbon advanced 

F 

F 

F 

F 

F 

F 

0 

0 


Greidanus 

F 

F 

F 

F 

F 

F 

0 

0 


Bourcier de Carbon elementary 

F 

F 

F 

F 

F 

0 

0 

0 


Melzer 

F 

F 

F 

F 

F 

0 

0 

0 

Assumes tire twist angle = swivel angle 

Moreland advanced 

F 

F 

0 

0 

F 

0 

0 

0 

Introduces time-constant term 

Moreland intermediate 

F 

F 

0 

0 

F 

0 

0 

0 

Implies extremely large value 

Moreland elementary 

00 

00 

00 

0 

F 

0 

0 

0 

Assumes laterally and torsionally rigid tire 

Temple elementary 

F 

00 

F 

F 

0 

0 

0 

0 

Assumes laterally rigid tire 

Maier 

F 

00 

0 

0 

0 

0 

0 

0 

Assumes laterally rigid tire 

Taylor 

00 

F 

00 

0 

00 

0 

0 

0 

Assumes torsionally rigid tire 

Kantrowitz 

F 

F 

F 

F 

F 

0 

0 

0 

J For trail not equal to zero both of these 

Wylie 

F 

F 

F 

F 

F 

0 

0 

0 

] theories can lead to erroneous conclusions 


-kche symbol F indicates a finite number. 
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Figure 1.- Geometric relations for a rolling elastic tire. (For the sake of clarity parts (b) 
and (c) of this figure are drawn to different scales. The positive X-axis is the direction 
of undisturbed motion.) H 
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(a) Assumed theoretical shape of tire equator distortion for a stationary- 

twisted tire, v = 0. 



(b) Actual shape of tire equator distortion for an untwisted tire at 
rest (solid lines) and just after starting to roll (dashed lines). 



(c) Actual shape of tire equator distortion for a rolling tire. 
Figure 2.- Tire equator distortion. 
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Figure 7.- Experimental stability boundary for a tire 29 cm in diameter 

(from ref. 17)* 
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Figure 8.- Qualitative comparison of the stability-boundary predictions 
for case I according to the summary and systematic approximation 
theories neglecting damper and gyroscopic moments (g = t = 0). 
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Figure 9 *- Comparison of theoretical and experimental predictions of the 
stability boundaries for the Von Schlippe-Dietrich test model of ref- 
erences 3 to 5 * (Hysteresis effects are neglected in these theoreti- rc 

cal curves.) ” 


NACA TN 3632 


Frequency, f, cps 



0 200 U00 600 800 1000 1200 


Velocity, v, cm/sec 

Figure 10.- Comparison of theoretical and experimental shimmy frequencies 
on the stability boundary for the Von Schlippe-Dietrich test model of 
references 3 to 5. (Hysteresis effects are neglected in these theo- 
retical curves.) 
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deal calculations of the 
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Von Schlippe-Dietrich theory (d-^ = 900 kg-cm/radian, d ? = 0.065 cm" 1 ) 

Approximation B (.r^ = 0.1, r) = 0.2 ) 

O Experimental data 
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Velocity, v, cm/sec 

Figure 13.- Effects of hysteresis on the shimmy frequency for the Von Schlippe- 

Dietrich test model of references 3 to 5* 
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tical trail, a . cm 



Inclination of swivel axis, k , deg 


Figure 14 . - Influence of swivel-axis inclination on the stability boundary. 
Tire approximately 12 cm in diameter; F z = 6.25 kg; v = 19 km/hr. 
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Figure 15.- Influence of shimmy amplitude on the shimmy frequency. 


Thebretical calculations 1 

Experimental data from reference 10 J Stability boundaries 


-Equation (lij.5) 

(hysteresis not considered) 


Equation ( lU7 ) 

^ (hysteresis considered) 


(Stable motion on right 
of stability boundaries) 


Velocity, v, cm/sec 


Figure 16.- Comparison of theoretical and experimental variations of criti- 
cal trail with rolling velocity, r = 3-5 cm; = 3-22 cm; N = 12.3 kg; 
= 5.9 kg-cm; r] = 0.1; = 0.2. 
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Dimensionless trail 


H 



Figure l8.» Variation of critical trail with rolling velocity according 

to Moreland's advanced theory. 
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Dimensionless damping constant 







Dimensionless damping constant, g/vmNi 



Dimensionless velocity, v Vm/Nr 


Figure 20.- Comparison of damping required for stability according to approximations C, Dl, 
D2, and D3 for a sample landing gear. L = 0.8r; h = a = 0.5rj c-j_ = c 2 = 0.25r; e = 0.3r> 

m-, = 0.35m; I Q = mr 2 j t = p = 0. 
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